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Abstract 

We calculate the 0{ag) heavy flavor contributions to the Wilson coefficients of the structure 
function F2{x, Q"^) and the massive operator matrix elements (OMEs) for the twist-2 oper- 
ators of unpolarized deeply inelastic scattering in the region massive Wilson 
coefficients are obtained as convolutions of massive OMEs and the known light flavor Wil- 
son coefficients. We also compute the massive OMEs which are needed to evaluate heavy 
flavor parton distributions in the variable flavor number scheme (VFNS) to 3-loop order. 
All contributions to the Wilson coefficients and operator matrix elements but the genuine 
constant terms at 0(q:^) of the OMEs are derived in terms of quantities, which are known 

(3) (3) 

for general values in the Mellin variable A^. For the operator matrix elements ^Qg, ^qg q 

the moments TV = 2 to 10, for ^gf ' to = 12, and for Afl^^"" , A^^^ 

to A^ = 14 are computed. These terms contribute to the light flavor -|— combinations. For 
the flavor non-singlet terms, we calculate as well the odd moments A^ = 1 to 13, corre- 
sponding to the light flavor — combinations. We also obtain the moments of the 3-loop 
anomalous dimensions, their color projections for the present processes respectively, in an 
independent calculation, which agree with the results given in the literature. 



^Present address: Institute de Fisica Corpuscular, CSIC-Universitat de Valencia, Apartado de Correros 22085, 
E-46071 Valencia, Spain. 



1 Introduction 



Deep-inelastic scattering processes of charged or neutral leptons off proton and deuteron targets, 
in the region of large enough values of the gauge boson virtuality = —q^ and hadronic mass 
= {q + pY, allow to measure the leading twist parton densities of the nucleon, the QCD- 
scale Aqcd and the strong coupling constant as{Q'^) = as{Q^)/{4:7i), to high precision. The 
precise value of Aqcd, a fundamental parameter of the Standard Model, is of central importance 
for the quantitative understanding of all strongly interacting processes. Moreover, the possible 
unification of the gauge forces [1] depends crucially on its value. Of similar importance is 
the detailed knowledge of the parton densities for all hadron-induced processes [2], notably for 
the interpretation of all scattering cross sections measured at the Tevatron and the LHC. For 
example, the process of Higgs-boson production at the LHC [3] depends on the gluon density 
and its accuracy is widely determined by this distribution. 

Let us consider the kinematic region in deeply inelastic scattering, where processes of higher 
twist can be safely disregarded and the hard scales and are large enough to allow the ap- 
plication of the light-cone expansion, saturated by the twist-2 contributions. The scattering pro- 
cesses are then described by structure functions -Fj(x, Q^), which decompose into non-perturbative 
mass/ess parton densities fj{x,ii'^) and perturbative coefficient functions Cl{x,Q'^ / ji^) by 

F,(x,Q2) U%)®fji^^f^") ■ (1-1) 

j=g,g,9 ^ 

The scale /j, denotes the factorization scale, which is arbitrary and cancels between the coefficient 
functions and parton distribution functions in the respective orders in perturbation theory. The 
symbol denotes the Mellin convolution 

[A®B]{x)= f dxi f dx26{x - xiX2)A{xi)B{x2) . (1.2) 
Jq Jo 

The Mellin transformation 

M[f(x)](N)= fdxx'^-'fix), (1.3) 
Jo 

if applied to (1.2), resolves the convolution into a product. 

Since we strictly consider twist-2 parton densities in the Bjorken limit, no transverse mo- 
mentum effects in the initial distributions will be allowed, which otherwise is related in the 
kinematic sense to higher twist operators. As is well known, the leading-twist approximation 
and the QCD improved parton model are equivalent descriptions for the dominant contributions 
to the deep-inelastic structure functions at sufficiently large scales Q^. The condition for the 
validity of the parton model [4] demands that 

— « 1 , (1.4) 

Tlife 

with Tint being the interaction time of the virtual gauge boson with a hadronic quantum- 
fluctuation, the life-time of which is given by rufc. The latter can be interpreted as a partonic 
state, provided (1.4) holds. Both times are measured in an infinite momentum frame and they 



2 



are given by 



1 APx 

Tint ~ — = jr^ r (1.5) 

1 _ 2P 



= (1.6) 



with P the large momentum of the hadron, the energy component of the virtual gauge boson 
in the infinite momentum frame, Ei the energy of the ith fluctuating parton, fc^ mj, Xi its 
transverse momentum, mass, and momentum fraction, E the total energy, and Mjv the nucleon 
mass. In the region of not too small values, nor values near the elastic region a; ~ 1, of the Bjorkcn 
variable the partonic description holds for massless partons. Evidently, iff Q^{1 — a;)^/mf 
no partonic description for a potential heavy quark distribution can be obtained. In the general 
kinematic region the parton densities in Eq. (1.1) are enforced to be massless and the heavy quark 
mass effects are contained in the Wilson coefficients Cf, which are perturbatively calculable. Due 
to this, one may identify the massless flavor contributions and separate the Wilson coefficients 
into a purely light part cf^^^"^^ and H^, which accounts for the heavy quark contributions, 

C/(.,f)^C/*'(.,f)+H;(.,f,^),*^c,6. (1.7, 

The question, under which circumstances one may introduce a heavy flavor parton density, will be 
discussed later. Both, the measurements of the heavy flavor part of the deep-inelastic structure 
functions, cf. [5], and numerical studies [6] based on the leading [7] and next-to-leading order 
(NLO) heavy flavor Wilson coefficients [8], show that the scaling violations of the light and the 
heavy contributions to (1.7) exhibit a different behaviour over a wide range of Q^. This is both 
due to the logarithmic contributions ln'^(Q^/m^) and power corrections oc {w? /Q'^Y, k > 1. 
Moreover, in the region of smaller values of x the heavy flavor contributions amount to 20- 
40%. Therefore, the precision measurement of the QCD parameter Aqcd [9] and the parton 
distribution functions in deeply inelastic scattering require the description of the light and heavy 
flavor contributions at the same accuracy. The separation (1.7) allows the definition of the 
light flavor contributions and the related heavy flavor contributions to Fi{x,Q'^) applying the 
factorization Eq. (1.1). 

The perturbative accuracy reached for fI^^'^^{x, Q^) is of 3-loop order [10-26], which requires 
to calculate the 3-loop heavy flavor Wilson coefficients as well. The NLO heavy flavor corrections 
in the complete kinematic range arc available only in semi-analytic form [8] due to the complexity 
of the contributing phase space integrals. ^ Heavy flavor corrections to different sum rules for 
deep-inelastic structure functions were calculated in [28]. An important part of the kinematic 
region is that of larger values of Q^. As has been shown in Ref. [29], the heavy ffavor Wilson 
coefficients h\2{x , Q'^ / fi^ , mf / fi^) can be calculated analytically at NLO for Q"^ /m'^ >10. ^ This 
is due to a factorization of the heavy quark Wilson coefficients into massive OMEs, Aj^, and 
massless Wilson coefficients, Cj'^^^^^ in case one heavy quark ffavor of mass m and Uf light ffavors 
are considered. This restriction to only one heavy quark ffavor is required beginning with the 
3-loop corrections and will be adopted in the following. In the present paper, we calculate the 



precise numerical implementation in Mellin space was given in [27]. 
•^In case of H]^(.;:, Q^/^^, m?//i^) this approximation is only valid for Q^/m^ > 800, [29]. The 3-loop corrections 



were calculated in Ref. [30]. 



3 



massive operator matrix elements Ajk contributing to the heavy flavor Wilson coefficients for 
the structure function F2{x,Q'^) in the region Q'^/m'^> 10 to 3-loop order for fixed moments of 
the Mellin variable N. In case of the flavor non-singlet (NS) contributions, we also present the 
odd moments of the — projection. We further calculate the operator matrix elements, which are 
required to define heavy quark densities in the VFNS [31]. Due to renormalization, higher order 
contributions in e to corrections of lower order in Og, cf. [29,31-35], and other renormalization 
terms, such as the anomalous dimensions and the expansion coefficients of the QCD /3-function 
and mass anomalous dimensions, contribute. For these reasons, the present calculation yields 
also the moments of the complete 2-loop anomalous dimensions and the terms oc Tp of the 
3-loop anomalous dimensions ■jij{N). In the pure singlet (PS) case, Iqq^^iN), and for ■jgg^N), 
these are the complete anomalous dimensions given in [18, 19], to which we agree. Since the 
present calculation is completely independent by method, formalism, and codes, it provides a 
check on the previous results. Except for the constant part of the unrenormalized heavy flavor 
operator matrix elements, we obtain the heavy quark Wilson coefficients in the asymptotic region 
for all values of the Mellin variable N. The analytic continuation of these expressions to complex 
values of N can be performed with the help of the representations in [36] and those given for 
the anomalous dimensions and massless Wilson coefficients in [18,19,26]. 

The paper is organized as follows. In Section 2, a brief outline of the basic formalism is 
given. The renormalization of the different massive operator matrix elements is described in 
Section 3. In Section 4, we present details on the unrenormalized and renormalized operator 
matrix elements. Technical details of the calculation and the main results are discussed in 
Section 5. Depending on the CPU time and storage size required, the moments up to = 10, 12, 
and 14 of the different operator matrix elements could be calculated. In Section 6, representations 
for heavy quark parton densities in the region /i^ ^ are given and Section 7 contains the 
conclusions. In the Appendices, we give a consistent set of Feynman rules for the composite 
operators up to 3-loop order, present the moments of the 3-loop anomalous dimensions, and of 
the constants part of the different 3-loop massive operator matrix elements. 

2 The Formalism 

The heavy quark contribution to the structure function F2{x, Q^) for one heavy flavor of mass 
m and Uf light flavors is given by, [31], 






(2.1) 
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with (S) the singlet contributions. Here, we denote the heavy flavor Wilson coeflicients H] by Lj, 
Hf respectively, depending on whether the photon couples to a hght (L) or the heavy (H) quark 
line. /fc(a;, yU^) and f-^[x,jjL^) denote the quark- and antiquark distribution functions, G{x,ii^) is 
the gluon distribution and 

S(x, ii') = J2 [M^, /^') + fki^, /^')] (2-2) 
k=i 

denotes the flavor singlet distribution, eg is the electric charge of the heavy quark. Due to the 
difference of quantities taken at n/ + 1 and n/ flavors, it is useful to adopt the following notation 
for a function /(n/), 

/(»/) = /("/ + !)-/(»/). (2.3) 
/(»/) - (2.4) 

and/(n/) = [/(n/)]. ^ As has been shown in Ref. [29], the heavy quark Wilson coefficients 

in deeply-inelastic scattering, H^, factorize in the region ^ m^, in which power corrections 
can be disregarded, into massive operator matrix elements A^f'^ and the light flavor Wilson 
coefficients C-/^', 

- AT'' ® c:^r ' (2-5) 

where i = 2,L specifies the structure function considered. 

The operator matrix elements A^f^ are the partonic expectation values 

A'if'' [n,^) = {l\Of''\l) , l-q.g, (2.6) 



with the local twist-2 operators given by 



A. 



OZ,,u...,,n = ^"-^S[^7m,Z}^2...Z}^„:^^]- trace terms, (2.7) 
0|;;ui,...,^„ = r"^S[^7^,D^2 . . . D^^ip] - trace terms , (2.8) 
Ol,,,,...,,„ = 2z-2SSp[F;^,D^, . . . - trace terms , (2.9) 

for the fermionic non-singlet, singlet, and gluonic case, [37]. Here, S denotes the symmetrization 
operator of the Lorentz indices /ii, . . . , Xa is the flavor matrix of SU{nf) with Uf light flavors, 
denotes the quark field, the gluon field-strength tensor, and the covariant derivative. 
Sp in (2.9) is the color-trace. The quarkonic operator matrix element can be represented by 

A\ = A'' + <'- (2-10) 

The different contributions to (2.5) were given in [31], Eqs. (2.31-2.35). To O(a^), the Wilson 
coefficients in MeUin space are : 



NS,(2)/ X . ANS,(2), 



^Later on, the symbol ' will also be used for the bare coupling dg, the mass m, and the bare OMEs, where 
(2.3) is not applied. 
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+ a 



NS,(1), 



VNS,(3) , 



.PS,(3) 



+ a 



+ a' 



A, 



PS, (3) 
Qq 

lPS,(2), 



(«/) + 



PS, (3), 



(^/+i)+4?,Q(^/)c^:i(^/+i) 



+ at 



+ at 



c'S(% + i) 



,NS,(1) 



(2) 



^PS,(2). 



.(2)/ 



(2.11) 
(2.12) 



(2.13) 



(2.14) 



+ (7g(n/ + l) 



(2.15) 



For brevity, we have dropped here part of the arguments of the Wilson coefficients and operator 
matrix elements by identifying = H^(A^, /x^/m^, n/), = Cf{N,Q^/n'^,nf ) and Aij = 

Aij{N, / /i'^ , n f) . These representations were verified in the LO and NLO case comparing with 
the results in [7,8] for Q"^ m^. 

The massive operator matrix elements are calculated keeping the external massless parton 
lines on-shell, while the heavy quark mass m sets the scale. The massless Wilson coefficients 
in (2.11-2.15) were calculated in Refs. [20,21,23-26]. 



Renormalization of the Massive Operator Matrix Ele- 
ments 



We perform the calculation of the massive operator matrix elements in D = 4 + £ dimensions 
and apply dimensional regularization. For each loop integral a factor 



exp 



ln(47r)) 



(3.1) 



with 7e the Euler-Mascheroni constant, is obtained which collects universal terms, and Ss ■= 1 in 
the MS-scheme. The following equation shows the pcrturbativc expansion of the unrenormalized 
OMEs, denoted by a double-hat, in the bare coupling constant a^ in Mellin space 



-^2^Z,/2,(0 

A.J I m- 



=1 



= /i^£,7v),(3.2) 
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with 



4.(m^ = ;.^s,A^)=^-l^a(f)(A^). (3.3) 

k=0 

Here, is the Melhn-parameter, (1.3), m the bare mass, and fi = fiR is the renormahzation 
scale. Also the factorization scale fip will be identified with n in the following. 

The factorization between the massive OMEs and the massless Wilson coefficients (2.5) re- 
quires the external legs of the operator matrix elements to be on-shell, 

= , (3.4) 

where p denotes the external momentum. Unlike in the massless case, where the scale of the 
OMEs is set by an off-shell momentum — < 0, in our framework the internal heavy quark mass 
sets the scale. In the former case, one observes a mixing of the physical OMEs with non-gauge 
invariant (NGI) operators, cf. [16,38,39], and contributions originating in the violation of the 
equations of motion (EOM). Terms of this kind do not contribute in the present case. 

The renormahzation of the massive OMEs is performed in four steps. First mass renormal- 
ization is carried out, for which we use the on-mass-shell scheme and later also compare to the 
results in the MS-scheme. Afterwards, charge renormahzation is performed in the MS-scheme. 
To maintain condition (3.4), which is of physical importance, we will, however, first introduce 
a MOM-scheme for the strong coupling constant and then perform a finite renormahzation 
changing to the MS-scheme. The former scheme is implied by keeping the external massless 
parton lines on shell. Note, that there are other, differing M M-schemes in the literature, cf. 
e.g. [40]. After mass and coupling constant renormahzation, the OMEs are denoted by a single 
hat, Aij. The ultraviolet singularities of the composite operators are canceled via the corre- 
sponding -factors and the UV-finite OMEs are denoted by a double tilde, Aj^. Finally, the 
coUinear divergences are removed via mass factorization. 



3.1 Mass Renormalization 

There are two main schemes to perform mass renormalization: i) the on-shell scheme and ii) 
the MS-scheme. We will apply the on-shell scheme in the following, defining the heavy quark 
mass as the pole mass, and compare to the MS-scheme later. The bare mass in (3.2) is replaced 
by the on-shell mass m through 



m 



ZmTn = m 



1 + 



^ /"m \ 



6mi + a 



S \ 9 



+ Oial 



(3.5) 



The constants in the above equation are 



5mi — Cf 



--4+(4+^C2U 



+ 5m?) + Sm?e , 

e 



1 



- (18Ci. - 22Ca + 87V(n/ + Nn)) + - -—Cf + —Ca - 147>(n/ + Nn) 



45 



91 



(3.6) 
(3.7) 



^Note that there is a misprint in the double-pole term of Eq. (28) in Ref. [34]. 
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199 51 

_ _ + 481n(2)C2 - 12(3 ) + Ca 



45 \ /69 



fine; c; 

■^ + |C2-241n(2)C2 + 6C3 



(5m^ (0) 
^ ^ h drrin 



(3.8) 



(3.9) 



with = (N^ - 1)/{2Nc),Ca = N^.Tp = 1/2 for SU{Nc) and iV^ = 3 in case of QCD. Ck 
denotes the Riemann (^-function. In (3.8), Uf denotes the number of hght flavors and Nh the 
number of heavy flavors, which we will set equal to one from now on. The pole terms were given 
in [41,42], and the constant term in [43,44], see also [45]. In Eqs. (3.7, 3.9), we have defined 
the expansion coefficients in s of the corresponding quantities. The following equation shows the 
general structure of the OMEs up to 0{dl) after mass renormalization 



4(^,.,iv) 



5ij + dsAij \^—,e,N 



(3)/^j2 . /m^\'^/2 d 5(2) /^2 



+5m.(^) 



dm ' '^ V ij? 



d^ 5(1) 



m 



dm 



m 



,e,N 



(3.10) 



3.2 Renormalization of the Coupling 



As the next step, we consider charge renormalization. We briefiy summarize first the main steps 
in the massless case in the MS-scheme. Afterwards, we extend the description to the massive 
case in the MOM-scheme which we use, before we transform back to the MS-scheme. 
The bare coupling constant dg is expressed by the renormalized coupling a^^ via 



(3.11) 



The coefficients in Eq. (3.11) are, [46-49] and [50,51], 



(3.12) 
(3.13) 



with 



4^ 



A(n/) 
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Cl-4(^^CA + CF^TFnf . 



(3.14) 
(3.15) 



The evolution equation for the renormahzed couphng constant is then given by 



(3.16) 



fc=0 



The factorization relation (2.5) strictly requires that the external massless particles are on 
shell. Massive loop corrections to the gluon- and ghost-propagators violate this condition, which 
has to be enforced subtracting the corresponding corrections. They can be uniquely absorbed 
into the strong couphng constant applying the background field method [52-54]. Here, Zg can be 
obtained by only considering the gluon propagator. After mass renormalization in the on-shell 
scheme via Eq. (3.5), we obtain for the heavy quark contributions to the gluon self-energy 

— I I exr> I > — I 

.2. 



+al(^y ^(^-'^TfCa-^TfCf)-^TfCa + 15TfCf 



86 



— TpC^ — —TfCf — -^C^TfCa — C^TfCf^ 



(3.17) 



Note, that although the 0(a5)-term in the above formula is an expression to all orders in £, the 
0(ag)-term and hence the formula in general only holds up to 0{e). We have used the Feynman 
rules of the background field formalism as given in Ref. [55]. In the following, we define 



exp 



The renormalization constant of the background field Za is related to Zg via 



^-2 



(3.18) 



(3.19) 



The light-flavor contributions to Z^, Za,i, can thus be determined by combining Eqs. (3.11) 
and (3.19). The heavy flavor part, Za,h, follows from the condition 



nH,BF(0, as, m^) + Za,h = , 



(3.20) 



which ensures that the on-shell gluon remains strictly massless. Thus we deflne the renormal- 
ization constant of the strong coupling with Uf light and one heavy flavor as 

Z^^^(s,nf + l,^i,m) = 



Za.i + Z 



1 1/2 



(3.21) 



and obtain 



Z^°^'(£,m,/.,n^ + l) 



l + ar7/^')[-(/3o(ri/)+/3o,Q/(£)) 



+a 



MOm2j-^2'j 



+ 'e{$y{(^^^Q + '(^W + ''^i%)] +0(.^ar^■^) , (3.22) 



with 



32 
y 

—-^TpCa — —TpCp — C2 f -TfC'a + ^fC'f ) ■ 



(3.23) 
(3.24) 

(3.25) 
(3.26) 

(3.27) 

+ 0{e') . 
(3.28) 

Since the MS-scheme is commonly used, we transform our results back from the MOM- 
description into the MS-scheme, in order to be able to compare to other analyzes. This is 
achieved by observing that the bare coupling does not change under this transformation and one 
thus obtains the condition 



The coefficients corresponding to Eq. (3.11) expressed in the MOM-scheme read 

MOM _ 



2/3o(n/) ^ ^Po,Qj^^^ 



e e 
+ 



(6,n, + l)ar(/.^) = Z-"-(£,m,/.,rz, + l)ar7/.^) . 



,MS/,,2 



M0M2/ 



,M0M/,,2\ 



The following relations hold : 



MOM 



(1) 

Q 



(3.29) 



af , (3.30) 



or. 



,MS 



MOM I MOM^ 



Kr - Kiirif + 1) I + ar^"^ I 5a,7- - 5a-(n, + 1) 



-25a^^!{nf + 1) 5a^^?^ - 5a^^!{nf + 1) ) + ©(a^^^) , 



(3.31) 



vice versa. Eq. (3.31) is valid to all orders in e. Here, a^^ = a^^{nf + 1). Applying the on-shell 
scheme for mass renormalization and the described MOM-scheme for the rcnormalization of the 
coupling, one obtains as general formula for mass and coupling constant renormalization up to 

0(«MOM^) 



M0m2 

s 



'^(2) /rn2\e/2 ^ j 



M0m3 



+a 



+51711 



;{3) x(l) /x(2) /'n-)2x 



;(!)■ 



(1) 



\ jj? / dm 



2x^/2 d 5(2) 



am 



(i 5m? /m^ 



m——A-^ + 
dm 



13 



2 V/i2 



m 



rf2 ^(1) 



(im 



,(3.32) 



where we have suppressed the dependence on m, e and in the arguments. 



^Here we corrected a typographical error in [34], Eq. (48). 
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3.3 Operator Renormalization 



The renormalization of the ultra-violet (UV) singularities of the composite operators is done 
introducing the corresponding Zjj-factors. We consider first the case of n/ massless flavors, 
cf. [56], 



^2 , . ^2 



<^(^,ar,^/,iv) = Z-''^'{a^\nf,e,N)A^^{-^,a'^\nf,e,N) (3.33) 

2 2 

A,(^,a?^,n/,7V) = Z-\a^,nf,e,N)Aij(^:^,a^,nf,e,N^ , i, j,l ^ q, g,{3M) 

for the non-singlet and singlet case, with p a space-like momentum. As mentioned before, we 
neglected all terms being associated to EOM and NGI parts, since they do not contribute in the 
renormalization of the massive on-shell operator matrix elements. The NS and PS contributions 
are separated via 

= + ^.V''' ' (3-35) 

A,, = + (3.36) 

The anomalous dimensions of the operators are then given by 



iX'^^f^N) = fi^lnZ^^^{a^\nf,e,N) , (3.37) 



7i,(af ,n/,7V) = Zr^\a^^nf, e, N) i,-^Zij{a^^nf,e, N) . (3.38) 
They can be expanded into a perturbative series as follows 

1=1 

where the PS contribution starts at 0{al). The anomalous dimensions are known for all N at 
LO, [57,58], and NLO, [10-16]. Fixed moments at NNLO have been calculated in Refs. [23-25] 
and the complete result has been obtained in Refs. [18, 19]. At the level of twist-2, they are 
connected to the splitting functions, [59], by a Mellin-transform ^ 

7lf (n/, N) = - r dzz^'-'P^fiuf, z) . (3.40) 
Jo 

In the following, we do not write the dependence on the Mellin-variable N for the OMEs, the 
operator Z-factors and the anomalous dimensions explicitly. Furthermore, we will suppress the 
dependence on e for unrenormalizcd quantities and Z-factors. From Eqs. (3.37, 3.38), one can 
determine the relation between the anomalous dimensions and the Z-factors order by order in 
perturbation theory. In the general case, one finds 



^Duo to our convention, Eqs. (3.37, 3.38), there is a relative factor of 2 between the anomalous dimensions 
considered in this work and Refs. [18, 19]. 



11 



(2) 



(3.41) 



The NS and PS Z-f actors are given by 

1 + a 



7NS/_MS „ \ 



_ (0),NS _2 I 1 .1 , 



<29 



(3.42) 



qMS^ I ^ (0) (0) 1 (1),PS I +^MS'M Ifl (0) (0) (0) 
1 2^2 '95 '99 ^ 9e '99 f ^ ""s 1 ^3 V '99 '99 '99 



1 



+ 1^(0)^(0) (0) ^ ^ .U, 
' /- /oq /oo /OQ ' A-^U /go /, 



2£ 



£n3 



(0)^X0) 



g '99 '99 '99 



99 '99 



+ 



1 /I 



(0)^,(1) 



r7ofl 7( 



99 '99 



1 



1 



,(0)^,(1),PS 



^ Qh9 '91 ' 2 3 



+ 



799 



(2),PS 



3£ 



(3.43) 



All quantities in Eqs. (3.41-3.43) refer to n/ light flavors and renormalize the massless off-shell 
OMEs given in Eqs. (3.33, 3.34). 

In the next step, we consider an additional heavy quark with mass m. We keep the external 
momentum artificially off-shell for the moment, in order to deal with the UV-singularities only. 
For the additional massive quark, one has to account for the prescription of the renormalization 
of the coupling constant we used in Eqs. (3.27, 3.28). The Z-factors including one massive 
quark are then obtained by taking Eqs. (3.41-3.43) at + 1 flavors and performing the scheme 
transformation given in (3.31). The emergence of 5a^'^^ in Zij is due to the finite mass effects 
and cancels singularities which emerge for real radiation and virtual processes at — > 0. Thus 
one obtains 



(0) 

MOMjij_ I M0m2 



MOM„,{0) 



+/5o7 



+ a 



M0M3 



+ 



1 (0) , 1/'_1^2 (0) _ ^^^(0) (0) _ 1 (0) (0) (0) 

g7ii lij ) ^ 7ii Polii lij qIu Ilk ikj 



(3.44) 



and 



(0),NS 

7-l,NS/„M0M „ 1 „MOm 799 , ^M0m2 



1/'_1^(1).NS _ ^ MOM (0),NS\ 
^\ 2 '11 '99 J 



^In Eq. (3.43) we corrected typographical errors contained in Eq. (34), [34]. 
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1 



NS , 1^(0),NS2\ , M0m3 1 /'_1^(2),NS _ c MOM (1),NS 



(0),NS 



1 



7-l,PS/ MOM ^ , -, \ 



a 



M0M2 



(0),NS'; 



(3.45) 



i/'_l^a).PS^ , If 1^(0) (0) 



,M0m3 



_^ MOM (1),PSA + Ifl (0) (1) + 1 (0) (1) ^ 1 (0) (1) 



(0)^(0) (0) _ :^ (0) (0) (0) 

99 '99 '99 g '99 '99 '99 



(3.46) 



The above equations are given for n/ + 1 flavors. One rederives the expressions for nj hght 
flavors by setting (n/ + 1) =: n/ and 5a^'^^ = da^^. As a next step, we split the OMEs into a 
part involving only light flavors and the heavy flavor part 



Ai(p^m^/.^ar^r^, + l) = ^,4^,ar,r^/)+ig(p^m^/.^ar°^n, + l) .(3.4?) 



In (3.47, 3.48), the light-flavor part depends on a^^, since the prescription adopted for coupling 
constant renormalization only applies to the massive part. Af^ denotes any massive OME we 
consider. The correct UV-renormalization prescription for the massive contribution is obtained 
by subtracting from Eq. (3.47) the terms applying to the light part only : 



~Q 



A;.(p^m^/x^ar°^n, + l) 



+Z-^(ar^n, + l,/.)4(^,af,n,) 



-Z^^(af,n;,/.)4(-|-,af,n, 



(3.48) 



where 



(3.49) 



k=l 



In the limit — 0, integrals without a scale vanish within dimensional regularization. Hence 
for the light-flavor OMEs only the term 5^ remains and one obtains after expanding in 

i:j(^.arr"/ + l) = ar"(ir(f)+Z-''''(n/+l./-)-ZS''<''("/) 
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The Z-factors at n/ + 1 flavors refer to Eqs. (3.'! 
to the massless case. 



-z^'''''(»/)+2.;'''"(% + i.A)4r(f) 

(3.50) 

-3.46), whereas those at n/ flavors correspond 



3.4 Mass Factorization 

Finally, we have to remove the coUinear singularities contained in A^j, which emerge in the limit 
= 0. They are absorbed into the parton distribution functions. As a generic renormalization 
formula, generahzing Eqs. (3.33, 3.34), one flnds 

A,- = Z-,'A,T^I . (3.51) 

The renormalized operator matrix elements are obtained by 

4(^,ar^n, + l) = i;?(^,ar^% + l)r-. (3.52) 

If all quarks were massless, the identity, [29], 

r,, = Z-/ . (3.53) 

would hold. However, due to the presence of a heavy quark Q, the transition functions r(nj) 
refer only to massless sub-graphs. Hence the F-factors contribute up to O(a^) only and do not 
involve the special scheme adopted for the renormalization of the coupling. Due to Eq. (3.53), 
they can be read off from Eqs. (3.41-3.43). 

The renormalized operator matrix elements are then given by: 



2 

MOM 



(f ) + + 1) - z^''«w]r,7''"(ny) 
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+ 



-1,(1), 



(3.54) 



From (3.54) it is obvious that the renormahzation of Af^ to 0{a^) requires the 1-loop terms up 

to 0{e'^) and the 2-loop terms up to 0{e), cf. [29,31,33-35]. Finally, we transform the coupling 
constant back to the MS-scheme by using Eq. (3.30). We do not give the explicit formula here, 
but present the individual renormalized OMEs after this transformation in the next Section as 
perturbative series in a^^, 

4(f.''f.»/ + i) = W(^.»/ + i)+"^'-*«"(^'"'+0 



+«f'4''='(^."/ + i). 



(3.55) 



4 General Structure of the Massive Operator Matrix El- 
ements 

In the following, we present the unrenormalized and renormalized massive operator matrix ele- 
ments for the specific flavor channels. The pole terms can all be expressed in terms of known 
renormahzation constants, which provides us with a strong check on our calculation. In par- 
ticular, we obtain the moments of the complete anomalous dimensions up to 0{al), as well as 
their Ti?-terms at 0{al). The moments of the 0(£'^)-terms of the unrenormalized OMEs at 

(3) (2) 

the 3-loop level, a^j , are a new result. Previously, the 0{e) terms at the 2-loop level, a.-^ , for 
general values of N were calculated by the present authors in Refs. [34, 35] . The pole terms and 
the 0{e^) terms, a\^\ at the 2-loop level have been calculated for the first time in Refs. [29,31]. 
They were confirmed in [33,35], as well as by the present calculation, in which they appear 
in the renormahzation of the respective moments of the 3-loop OMEs. In order to keep up 
with the notation used in [29,31], we define the 2-loop terms al^\ a^f after performing mass 
renormahzation in the on-shell scheme. This we do not apply for the 3-loop terms. We choose 
to calculate one-particle reducible diagrams and therefore have to include external self-energies 
containing massive quarks into our calculation. Before presenting the operator matrix elements 
up to three loops, we first summarize the necessary self-energy contributions. 

4.1 Self— energy contributions 

The gluon and quark self-energy contributions due to heavy quark lines are given by 

tlf,{p\m\ ix", as) = iS'^' [-g^^p" + p^p,] Ii{p\m\ ii\ a,) , (4.1) 

with 

oo 

^(p^m^/.^a,) = J]a,^^(^)(p^m^/i^). (4.2) 

fe=i 
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and 

ilijip^^m^^l/^as) = i Sij :^E{p^,m'^,iJ,^,ds) , (4.3) 

where 

oo 

^:(p^m^/.^a,) = ^a,^EW(/,m^/.2). (4.4) 

k=2 

Note, that the quark sclf-cncrgy contribTitions start at 2-loop order. These self-energies are eas- 
ily calculated using MATAD, [60], cf. Section 5. The expansion coefficients forp^ = of Eq. (4.2, 
4.4) are needed for the calculation of the gluonic and quarkonic OMEs. The contributions to 
the gluon vacuum polarization for general gauge parameter ^ are 

"'■'(»•$) --^(f)"(-|-(|f(D'))- («) 



+H I ^ + - I) + C - 3G I M + O(s') (4.6) 



■ I ^ {Ca - 6Cf) n,TF + ^ (sSCj - 48Cp) Tp + ^| (-781 + 63^) 
+^C.,Cf I + j|^(Ca(-101 - I8C2) - 62CF)n,Tp 



SO 



3203\ 



-K-fc.-^)}.^.{c.(-f|..f..^i|) 



42799 
3888 



95 274 

y^^ + "9" 



+ 0(£) , 



(4.7) 
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and for the quark self-energy, 



^ 2 

m 



89 C2 
72 2 



2 X 3e/2 



- 0{e') (4.8) 
Ca{1 - a + ^^(^/ + 2) - C^(f + 4e) - ^C^} 

r8 604^ rl7 5 1879 7 3 407 1 
+Ci.{-8C3-C2-^}j+0(£), (4.9) 



see also [40,61]. In Eq. (4.7) the constant 



B4 = -4C2ln'(2) + -ln"(2) 
K -1.762800093... . 



-fa + i6m(i) 



(4.10) 



appears. 



4 2 A^'=> 



The lowest NS-contribution is of O(a^), 



A 



NS 



a 



(4.11) 



The expansion coefficients are obtained in the MOM-scheme from the bare quantities, using 
Eqs. (3.32, 3.54). After operator renormalization and mass factorization, the OMEs are given 

by 



A 



.(2),NS,M0M 

qq,Q 

{3),NS,M0M 
qq,Q 



= A 
= A 



(2) ,NS,M0M 

qq,Q 

(3) ,NS,M0M 

qq,Q 



+ Z-/'(^)''^^(n, + l)-Z-W^(n,), 
+ Z-/'(^)'-^(n, + l)-Z-/'(3)-NS(^^) 



(4.12) 



(2),NS,M0M 



+ 



^(2)^NS,M0M ^ ^-^i,(2),NS(^^ ^ _ z '^^^ (u f)] r-'W(n,) . (4.13) 



Prom (3.32, 3.54, 4.12, 4.13), one predicts the pole terms of the unrenormahzed OME. At second 
and third order they read 



j(2),NS 
A 

^qq,Q 



/m 



n (0) ^(1),NS 

£2 ^ 2e 



+ a(2).NS , ^(2),NS 



(4.14) 
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X(3),NS 
A 



/m \ 



3£3 



(/?o + 2/?o,q) + 



(1),NS 



0/3 (0) 



1 / -(2),NS 



/3o + Pq,q 



(3).NS 



99:' 



(4.15) 



Note, that we have already used the general structure of the unrenormahzed lower order OME 
in the evaluation of the 0(ag) term, as we will always do in the following. Using Eqs. (4.12, 
4.13, 3.32), one can renormalizc the above expressions. In addition, we finally transform back 
to the MS-scheme using Eq. (3.30). Thus one obtains the renormalized expansion coefficients 
of Eq. (4.11) 



A 



(2),NS,MS 
qq,Q 



A 



(3),NS,MS 

qq,Q 



Q (0) 

(^G^QTqq fm 
4 X/i 



^(l),NS 



/3 (0) 
{2),NS _ PO,Qlqq 

qq,Q A 



ln^(f ) + ||27r^/?o,Q - 27r'(/5o + /3o,c, 



Po + 2/3o,, 



(2) 7qqPoPo,QC3 



MS 

Q 



(3).NS 
^99,<3 



(4.17) 



Note that in the NS-case, one is generically provided with even and odd moments due to a 
Ward-identity relating the results in the polarized and unpolarized case. The former refer to the 
anomalous dimensions 7^,^'^ and the latter to Jqq'" as given in Eqs. (3.5, 3.7) and Eqs. (3.6, 
3.8) in Ref. [18]. The relations above also apply to other twist-2 non-singlet massive OMEs, as 
to transversity, for which the 2- and 3-loop heavy flavor corrections are given in [62]. 

^'^ ^Qq ^qq,Q 

There are two different PS-contributions. The term referring to the case in which the operator 
couples to a heavy quark, starts at 0(^5), whereas the term in which it couples to an 

internal light quark line, A^^q, emerges for the flrst time at O(a^), 



4 PS 
^Qq 
4 PS 
"^99,Q 



Qq 
4\ 



(4.18) 
(4.19) 



Separating these contributions is not straightforward, since the generic renormalization formula 
for operator renormalization and mass factorization, Eq. (3.54), applies to the sum of these terms 
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only. At 0{al), this problem does not occur and renormalization proceeds in the MOM-scheme 
via 



^^2),PS,M0M ^ ^g),PS,MOM^^_i,(2),PS(^^^l)_^-^l,(2),PS 



qq 



qq 



[n 



f) 



+ 



Thus the unrenormalized expression is given by 

;(2),PS 



A 



Qq 



- 2 . / -(0) (0) -^(1),PS 

m'y iWiW ifW (2),ps , _(2),ps^ 



(4.20) 



(4.21) 



The renormalized result in the MS-scheme reads 



(2),PS,MS 



— (-j+^ln(-j+a^, +^(2- (4.22) 



^Qq - 8 ^^^V//2;- 2 

The corresponding renormalization relation at third order is given by 

^g),PS,MOM ^ ^(2)^PS,M0M ^ ^g,PS,MOM ^ ^(3)^PS,M0M ^ ^^-^1,(3),PS (^^ ^ _ Z-.(3),PS 



-^."/'^'Hn/) 

_ v-l,(2),PS| 

qq ' 



r:i'(^)(n/) + 



39 



i(2),M0M 



(4.23) 



Taking into account the kinematic and UV-structure of the contributing Feynman diagrams, the 
two contributions can be separated. For the bare quantities we obtain 



X(3),PS 



^2 3^/2 



47i?'^^ 



(0) (1) (0) 

7gg Iqg _|_ 7gg 



2^(1) _ ^(1) 

'59 '59 



1 / -(2), PS 



+TS'4?,«-TS'<'gj-''(/'o+/'«.«)«?r- 



Iqg Igq ^2 
16 



e \ 3 

'^99 



i(2),PS 
Iqq 



99 



(4.24) 



;(3),PS 



99,< 



/m2^^3£/2 

, /~(2),PS 

e \ 3 



^799 739 PO; 



Q 



+ 7™ a 

' Ian ^ 



3^3 3^2 1 

:,(0).,(0) 



(2) 795 759 /5o,qC2 



95 59,<3 



+ 



J3),PS 



(4.25) 
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The renormalized terms are given in the MS-scheme by 



,(3), PS, MS 



799 Igq J (0) 
'99 



48 



1 
8 



+^(0) ("^(1) 

' '99 \ '99 '99 y 



^(0)^(1) 
'99 '99 



^(2), PS 
99 



-32a 



(2),PS 
Q9 



+ /?n + 8-7(0)0^^^ - 8-7(0) r;^^^ _ ^(0) (0) A /^^(o) _ 

l,PO -t- PO,Qj + ^Iqg Ugq^Q 07_w "Q(, 7gg Iqq S,2 ^ Tflfl 7, 



99 (39 '99 '99 



'99 '99 



^99,(3 



+6/3o + 8/3o,q) I In(^) + 4(/3o + /3o,Q)agf ' + 7S^agl - "fiS^ 

(0)-(0)> -(0) (1)> 

^799^1^ ^.0) _ .0) + A + 7,9^ _ ^^(1)^(0)^(0) + SmfHr^ 



, or (-1) (2), PS , (3), PS 



A 



(3), PS, MS 
qq,Q 



= Uf 



"Qq 
(0) ^ (0) ^ 

r9qrq9pO,Q 



(4.26) 



12 



^99 S9,Q 



7ot 7J?/3o,qC3 799^'''^/3o,qC2 I , (3),PS 



s«-7S'/3o,«c4)i.<^ 

X,(1).PS 



12 



+ a 



qq,Q 



(4.27) 



4.4 Agg and Agg^Q 

The OME is the most complex expression. As in the PS-case, there are two different 
contril:)utions, depending on whether the operator couples to a light quark line, denoted by 
Aqg^Q, or to a heavy quark hne, given by Aqg, 



Aqg^Q = alAf^^Q + 0{at) . 



In the MOM-scheme the 1- and 2-loop contributions obey the following relations 



(4.28) 
(4.29) 

(4.30) 



.(2),M0M _ 1(2), MOM 7-1,(2)/ 7~W) \ , y-hW ,xi(l),MOM 

^Q9 - ^Q9 + ^99 y^f + ^)- ^qg y^f) + ^qg l^/ + ^)^g9,Q 



;(1),M0M 



+ 



igr°^ + Z-/.«(n, + l) 



-Zqg''^'\nf)y-l'^'\nf) 
The unrenormalized terms are given by 

(f)"fexp(E|(£)-), 



(4.31) 



^Qg 



(0) 



(4.32) 



1=2 
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^(2) 



/my 



V 



-(0) 

2e2 \ '99 



Iqq ^ 



2/?o + 4/?o,q) 



7gg -26mi (2) _ c (0)^(0) 



2e 



(1)^(0) _ 7gg Po,q(,2 
12 . 



(4.33) 



Note that we have aheady made the one-particle reducible contributions to Eq. (4.33) explicit, 
which are given by the 1-Ioop contribution multiplied by the 1-loop term of the gluon-self 
energy, cf. Eq. (4.5). Furthermore, Eq. (4.33) already contains terms which result from 
mass renormalization in the 0{e^) and 0{e) expressions. At this stage of the renormalization 
procedure, they should not be present, however, we have included them here in order to have 
the same notation as in Refs. [29, 31] at the 2-loop level. The renormalized terms then become 
in the MS-scheme 



.(1),MS 



.(2),MS 



'(0) 

2 \^ 



(4.34) 



.(0) 

Tqg 
8 



-(0) A 

+«S + ^(tS'-tS'+2/J„). 

The generic renormalization relation at the 3-loop level is given by 

4r°" + ^fj:^°" = 4'r°" + ^s'i""" + z«'''Hnf + 1) - 



•(1) 

1 \ Tqg 1 fm'-\ 

'"(^)+i^'n^) 



(4.36) 



MOM 



+ Z-'(^)(n, + l)igf + + ^-/■(^)(n, + 1) - Z-/'«(%)lr-'(^H-/) 



+ 



- ^.-/'^^^'^^(-/)]r,-^'«(n,) + [<f '^^^ + ^-/■^'^■''(-/ + 1) 
-^,-/'^'^'''(^/)]r,-/'^^H%)- (4.36) 

Similar to the PS-case, the different contributions can be separated and one obtains the following 
unrenormalized results 



i(3) 



m 



2x 3e/2 



■-(0) 

Tqg 



7, 



(°) - 27(5) - 6/3o - 8/3o,Q + 8/3o' 



99 



+28/5o,q/5o + 24/3o^,Q + 7(5) 7^^^^ + 6/5o + 14/5o,q 



.(0) 



+ 



1 

6^2 



^(1) 

igg 



N?-27S^-8/5o 



+ ^(0) 



-10/3o,q] 



7S)''^^{1 - 2n,} + 7«'^^ + 7j?'^^ + 27^) - 7^^ " 2^ - 2/?,,, 
- 7j? + 3/3o + 5/3o,Q 



'''sfl '99 



+ £ 3 3 ^^^^ 



(2) (2),PS 
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(2) 

Qg 



7, 



(0) _ ^(0) _ 



1. 



99 



4/9o-4/?o,Q 



+ 



16 



(0) 



^(0) 
199 



6/?o + 2/3o,q} 



'99 



[-27; 



(3) 



(4.37) 



6^3 



7^£^ + 2/?o.c,7;^-7? + 2/?o 



^ r2 I 6 



2^(1) 



+7, 



«'^'-27jr' + 4/3i.Q| + 



91? 



+ i 2^ + ^(0) 



16 



7SHS^+2^o,g{ 



^(0) _ ^(0) 
199 '11 



+ 2/3o}] 



+ 



(2),NS 



(4.38) 



The renormalized expressions are 
A 



40) 

,(3), MS _ lqg_ 
^Q9 
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I (n/ + 1)7S^7S^ + 7S' (7^? - 27j? + 6/?o + 14/3o,q) + if^ {if^ - Qpo 



8/?o,q) + 8/3o' + 28/3o,Q/3o + 24/?^^ | In^ (^^) + i W (^7^^ 

^ 7S^ - 7^^ + (1 - n^nr' + + 7^^''' - 2A 

i(2) (2) 

-%^ + ^(7j?-7i?-4/3o-4/5o,Q) 



-6/5o, 



<3 ) 7q9 V 1 99 1 99 
2 



i,<3 



In^ 



v(2) 

;<39 



-(0) 
+ 7^f,(2) 



,(2),PS^ 



- 'V«S" -) + ^ (-('V + 1)7S'7S' + 7™ [27<? - 7S' - 6/3c 



2 



4/?o[2/3o + 3/3o,q]+7 



(0) 

99 



'99 ~ 



99 '99 

2. 



6/3o + 4/?o,c,]) ln(^)+ag(7S) 



-(0)/- 

-7j? + 4/?o + 4/?o,g) + 7J? (n/agj'" - ^L) + ^ ((-/ + H'J'f^ 



(0) 

99 



7^? - H'J + 6/3o - 2/3o,Q 



~ 199 



~{l)n /- -(0)> 
/99P0,qC2 799 C2 / (1) 
^ ^ ~ t \ ] nn li 



^ '99 



7j? - 6/3o 



+ 8/3o'-4/3o/3o,q-24/3o',q) 



(-1) 



16 

'Qg + '99 



16agl + 7(°) 



-24SmP - 8Smf^ - (2/^0 - 9(2/^0,^]) +^(27^^ 
-5m?^fj) + 5m?i'J (7^? - 7^? - 2/3o - 4/3o,q) + Smf^^'J + ag . (4.39) 
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4.5 Agq^Q 

The g'g-contributions start at O(a^), 



A 



,2/1(2) 



(4.41) 



The renormahzation formulae in the MOM-scheme read 



A 



(2), MOM 

9q,Q 



^Sr' + ^.V'^^Hn/ + 1) - ^-'(^n-/) 



+ A 



|(1),M0M 



(1) 



(4.42) 



j(3),M0M 



A^^^r^' + Z-.(3)(n, + 1) - Z-'(3)(n,) + Z-'(^)(n, + l)i2!?°" 



9q,Q 



+ 



(1),M0M 



j(l),MOM 



while the unrenormalized expressions are 



(4.43) 



i(2) 



i(3) 



7^ 



2/9o,Q 



■,(1) 



(0) , 7fl£_ (2) (2) 



(0) 



^2^ 3^/2 j ^^U, 
2 



7^ 



3£3 



2£ 

^(0)^(0) , 



91,Q 



''9q,Q'' 



(4.44) 



7, 



(°)-7S) + 10/3o + 24/3o,Q /9o,Q 



+ - 1 7i;^/3o, 



7i? - 7j? - 4/5o - 6/3o,Q 



+ 



(0) 

Tgq 



7 



'jr^ + 7S^''^^-7S^ + 2/?i,Q 



-45mS-^)/5o,Q7f°^ 



99 



1 / -(2) 
1 / 799 , i2) 



7S^ - 7j? - 6/?o,Q - 4/3o 



^ '99 



(2),NS 



"99,<9 + "<39 



(2), PS 



^99,Q 



, (0)^(1) , 799 C2 

+ 759 ^i,Q + 8 
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(4.45) 



The contributions to the renormahzed operator matrix element are given by 

#^^1^) + f ^K^) - , (4.46) 

^0 - 6/3o.q) + 7^? {^r' + 7^^'" - 7S^ + | In^ (^) 

,(2) 



.(2),MS 
.(3), MS 



+7S^(7S^-7j?-4/3o 



+'3!;i)+7S'C.(7<?7™ 



+«S«(7S'-7S'+''fl 

24 



J9q C3 (^^(0)^(0) , 



7j? - i'J + 12/3o,Q + lOPo 
^o + 6/3o,,)+7j?(4S.-^S^"-4r 



7j?-7S^ + 10/?o]/5o,Q 

gq + • 



799 Pl,Q 



(4.47) 



4.6 

The gfgf-contributions start at 0{as), 

A9,Q = ^sA% + a',AflQ + alAflQ + 0(at). 
The corresponding renormahzation formulae read in the MOM-scheme 

,M0M 

.Q 

[(2), MOM 



.(1),M0M 
.(2), MOM 



4»',r" + Z«'''"(n/+l)-Z„'-<"("/). 



99 V"./ ' / gg 

1,(1) 



(4.48) 



(4.49) 



(1),M0M 



+Z;/^^(n/ + l)A 



{1),M0M 

Qg 



+ 



4;!?°" + ^;9'''\-f + 1) - ^9-/'^^Hn,)lr-^(^)(n,) , 



(4.50) 



/.(3),M0M 
^99,Q 



i(3),M0M 
^99,Q 



99 



+ 



(1),M0M 

9g,Q 



A 



(2), MOM 

gg,Q 



+z--'<')(«, + i)<r1ri''">("/) 



(1),M0M 



+ 



4?,?°" + ^9-/'^'Hn/ + 1) - r-/'(^)(n,) . 



(4.51) 



24 



The general structure of the unrenormahzed 1-loop result is then given by 



99,* 



(0) 



^2 



An explicit calculation reveals 



4. ^ (^)"X-^)exp(Ef©')- 

^ i=2 



(4.52) 



(4.53) 



Using Eq. (4.53), the 2-loop term is given by 

X(2) 



99,Q 



1^' 



^{7S^7£^+2/3o,g(7S^ + 2/3o + 4/3o,Q)} 



(-1), 



+ 



2e 



,(2) ^OX^(0)« _^/:^2 . ^(2) , oX^W- .^o'qCs 



6 



(4.54) 



For Eq. (4.54) the same as for Eq. (4.33) holds. We have already included one-particle reducible 
contributions and terms stemming from mass renormalization in order to refer to the notation 
of Refs. [29,31]. The 3-loop contribution becomes 



A 



^(3) 



99,' 



ye/2 




(0) . (0) ^ 




£3 \^ 


6 I 



2/9o + 7(3o,Q 



6 



+ 



2/?2 + 7/3o,Q/3o + 6/3o',Q 

4/3o + 7/3o,Q 



(0) (1) - (1) 
73<Z759 739 



2/3o. 



Q 



-6/3o,q/?o 



(O).(O)^ 

16 



1 /-(2) 
1 / 799 



o (0)^ 



+ 



+ 



- 2(2/3o + 3/3o,Q)afiQ - nn^Q + + (^1,13% 



(0)^(2) 



.(0)„(2) ^ (0) 
'99 



7£)-7ff + 2(2n/ + l)/3o,Q + 6/3o 



+ 



Po,qC2 



7SH2/3o-/3o,q} + 4/3o' 



+5mf^ [-7;S(? - 27S)/3o,Q - 4/3o,q/?o - 



+ 25mt^(3o,Q 1 + afl^ 



. (4.55) 



The renormalized results are 



(4.56) 



l(2),MS 
99,Q 



7S^ 



8 



2/3o,g(7S^+2/5o)+7SHS^ 



+ a 



(2) 

'99,Q ' 



(4.57) 
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.(3),MS 
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Igq Iqg \ iqq ^gg 



6/?o - Anf(5o,Q - 10/?o,q) - ^(-ff^ [2/?o + 7/3o,q 
+4/3o^ + 14/3o,Q/3o + IWIq^O-Q^ In-^(^) + ^jTj? (tS^ + (1 - rz/)7^;0 
+7S^7S^ + 47S^/3o,Q - 47^) [/3o + Wo,q\ + 27^? + 4[/3i + (5,,q][3o,q | In^ (^) 

+7SHi?C2 (7^? - 7^? + 6/3o + 4ny/5o,Q + 6/3o,q) 

+4/3o,qC2(7S^ + 2/?o) (2/3o + Wo,q) | In(^) + 2(2/3o + S/^o.q)^?,^ 

(0)^(0)> ^ 

-6/?o) + ([/?o,Q - 2PoH'J + 2[(3o + 6/?o,q]/?o,q - ^P^) 

(tS^ + + ^ (7« - 27« - 2A - 2A,,) + ^ (e 



(2) 



+6m['^ (7j?7^? + 2/?o,Q7S^ + ^/^o.q/^o + 8/3o',q) - 25mf /3o,q + 4^,^ . (4.58) 



5 The Calculation of the Operator Matrix Elements 

In this chapter, we describe the computation of the 3-loop corrections to the massive operator 
matrix elements in detail. Typical Feynman diagrams contributing to the different channels 
are shown in Figure 1, where ® denotes the corresponding composite operator insertions, (2.7- 
2.9). The generation of these diagrams with the FORTRAN-based program QGRAF, cf. [63], is 
described in Section 5.1 along with the subsequent steps to prepare the input for the FORM- 
based program MATAD, [60]. The latter allows the calculation of massive tadpole integrals in D 
dimensions up to three loops and rehes on the MINCER algorithm, [64,65]. The use of MATAD and 
the projection onto fixed moments are explained in Section 5.2. Finally, we present our results 
for the fixed moments of the 3-loop OMEs and the fermionic contributions to the anomalous 
dimensions in Section 5.3. The calculation is mainly performed by using FORM programs, [66], 
while in a few cases codes have also been written in MAPLE. 



5.1 Generation of Diagrams 

QGRAF is a quite general program to generate Feynman diagrams and allows to specify vari- 
ous kinds of particles and interactions. Our main issue is to generate diagrams which contain 
composite operator insertions, cf. (2.7-2.9) and appendix 8.1, as special vertices. 

To give an example, let us consider the contributions to ^q^- Within the light-cone expan- 
sion, [67], this term derives from the Born diagrams squared of the photon-gluon fusion process 



26 




(NS) (PSh) (PSi) (qgH) (qgi) (gq) (gg) ghost 

Figure 1: Examples for 3-loop diagrams contributing to the massive operator matrix elements: NS - 
non-singlet, PSh,i - pure-singlet, singlet qgH,i, gq, gg and ghost contributions. Here the coupling of 
the gauge boson to a heavy or light fermion line is labeled by H and I, respectively. Thick lines: heavy 
quarks, curly lines: gluons, full lines: quarks, dashed lines: ghosts. 



shown in Figure 2. After expanding these diagrams with respect to the virtuahty of the photon, 
the mass effects are then given by the diagrams in Figure 3. These are obtained by contracting 
the hues between the external photons. Thus, one may think of the operator insertion as being 
coupled to two external particles, an incoming and an outgoing one, which carry the same mo- 
mentum. Therefore, one defines in the model file of QGRAF vertices which resemble the operator 
insertions in this manner, using a scalar field 0, which shall not propagate in order to ensure that 
there is only one of these vertices for each diagram. For the quarkonic operators, one defines the 
vertices 

(f) + (j) + q + q + n g , 0<n<3, (5.1) 

which is illustrated in Figure 4. The same procedure can be used for the purely gluonic interac- 
tions and one defines in this case 

(j) + (j) + n g , <n<A . (5.2) 

The number of diagrams we obtain contributing to each OME is shown in Table 1. The next 
step consists in rewriting the output provided by QGRAF in such a way, that the Feynman rules 
given in Appendix 8.1 can be inserted. Thus, one has to introduce Lorentz and color indices 
and align the fermion lines. Additionally, the integration momenta have to be written in such a 
way that MATAD can handle them. For the latter step, all information on the types of particles. 




Figure 2: Diagrams contributing to Aq^ via the optical theorem. Wavy lines denote photons; for the 
other lines, see Figure 1. 
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(1) 



Figure 3: Diagrams contributing to Aq^ 




<l> ■P2 




Figure 4: Generation of the operator insertion. 



the operator insertion and the external momentum are irrelevant, leading to only two basic 
topologies to be considered at the 2-loop level, which are shown in Figure 5. Note, that in the 
case at hand the topology on the right-hand side of Figure 5 always yields zero after integration. 
At the 3-loop level, the master topology is given in Figure 6. From this topology, five types of 
diagrams are derived by shrinking various lines. These diagrams are shown in Figure 7. After 
assigning the loop momenta, the Feynman rules are inserted. The computation of the Green's 
functions, which are associated to the respective operator matrix elements, still contain trace 
terms and require the symmetrization of the Lorentz indices. It is convenient to project these 
terms out by multiplying with an external source 



J. 



N 



(5.3) 



with being a light-like vector, A^ = 0. Additionally, one has to amputate the external field. ^ 
The Green's functions in momentum space corresponding to the local operators defined in Eqs. 
(2.8, 2.9) between gluonic states are then given by 



JNiA^ip) I o^-'^-(o) I Kip)) , 

Jn{a;{p) I o^-^no) I aI{p))q 

^'ip)Gt'Q,,u^''{p) = Jn{A;{p) I O^^-'^-(O) I Alip)) , 



(5.4) 
(5.5) 
(5.6) 



cf . [29] , with an external gluon field with color index a, Lorentz index /i, momentum p, and 
e^{p) the gluon polarization vector. In the flavor non-singlet case, Eq. (2.7), only one term 
contributes 



u 



{p,s)Gf'''XMp,s) 



JN{^^{P) I O,^--^-(0) I ^^{p))c 



(5.7) 



^Note that we choose to renormalize the mass and the couphng multipUcatively and thus have to include 
self-energy insertions containing a massive line on external legs. 
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Term 


# 


Term 


# 


Term 


# 


Term 


# 


.(3) 

^Qa 


1358 


.(3) 

^qg,Q 


140 


4(3),PS 


125 


4 (3), PS 

^qq,Q 


8 


4(3),NS 


129 


.(3) 
^9q,Q 


89 


4(3) 
^99,Q 


886 







Table 1: Number of diagrams contributing to the 3-loop heavy OMEs. 




Pi 



P2 



Figure 5: 2-Loop topologies, indicating labeling of momenta. 



with u{p, s), u{p, s) being the bi-spinors of the external quark and anti-quark, respectively. The 
remaining singlet and pure-singlet Green's functions with an external quark are given by, [29], 



u{p, s)G'^'^u{p, s) 
Hp, s)GI'qu{p, s) 
Hp, s)GI'qu{p, s) 



J^{^,ip) I o^-^-(o) I ^\p)), 
JN(^^{p) I oi;^-^-{o) I ^^{p))q 
Jn{Mp) I o^^-^-(o) I ^^{p))q 



(5.8) 
(5.9) 
(5.10) 



Note, that in the quarkonic case the fields \E', \1' with color indices i,j stand for the external 
light quarks only. The above tensors have the general form, cf. [29,56], 



Gab 



A 



m 



e,N]6''\A-p) 



- 9fiu + 



A 



1 = Q,g,q 



r,tj 



A,(^,e,iVj5^^(A-p)^^-^A, l = Q,9,q, r = S, NS, PS 



(5.11) 
(5.12) 



Here,yljj are the massive OMEs which we will calculate. In order to simplify this calculation, 
it is useful to define projection operators, which, applied to the Green's function, yield the 
corresponding OME. In the gluonic case, one defines 



ab;g l,fj,u 

p(2),fj.v^ab 
ab;g ^l,fiu 



Jab 



Snh 



ID 
1 



^ -(A-p)-^Gf 



(A-p) 



-N 



'9 



^ p>^A^+p^ 
+ A^ 



(5.13) 
(5.14) 
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Figure 6: Master 3-loop topology for MATAD, indicating labeling of momenta. 




In the quarkonic case, there is only one projector 

P,,,,Gr = ^(A.p)-^lTr[|^Gr]. (5.15) 
The unrenormahzed OMEs are given by 

= Piif'-Gj,, (5.16) 

K{^.e.N) = P,,,^ . (5.17) 

These projections yield the advantage that one does not have to resort to complicated tensorial 
reductions. In perturbation theory, the expressions (5.16, 5.17) can then be evaluated order by 
order in the coupling constant by applying the Feynman rules given in Appendix 8.1. While 
the projector (5.13) includes unphysical transverse gluon states, which have to be compensated 
adding the corresponding ghost- diagrams, (5.14) projects onto the physical states. 

To calculate the color factor of each diagram, we use the program provided in Ref . [68] . Up to 
this point, all operations have been performed for general values of Mellin N and the dimensional 
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parameter e. The integrals do not contain any Lorentz or color indices anymore. In order to use 
MATAD, one now has to assign to A?" a specific value. Additionally, the unphysical momentum A 
has to be eliminated by applying a suitable projector, which we define in the following section. 

5.2 Calculation of Fixed 3-Loop Moments Using MATAD 

We consider integrals of the type 

7Kp,m,ni,...,n,) ^ 1 | ^(A.g^r ■..(A.5,r^7(A;i, m) . (5.18) 

Here p denotes the external momentum, = 0, m is the heavy quark mass, and A is a light-like 
vector, A^ = 0. The momenta qi are given by any linear combination of the loop momenta ki 
and external momentum p. The exponents rii are integers or possibly sums of integers, see the 
Feynman rules in Appendix 8.1. Their sum is given by 

j 

^m^N . (5.19) 

i=l 

The function / in Eq. (5.18) contains propagators, of which at least one is massive, dot-products 
of its arguments and powers of m. If one sets = 0, (5.18) becomes 

7,(p,m,0,...,0)=7,(H = / 1^/(^1' -'^^H ■ (5.20) 

From = it follows, that the result can not depend on p anymore. The above integral is a 
massive tadpole integral and thus of the type MATAD can process. Additionally, MATAD can 
calculate the integral up to a given order as a power series in p^ /w?. Let us return to the general 
integral given in Eq. (5.18). One notes, that for fixed moments of N, each integral of this type 
splits up into one or more integrals of the same type with rij being just integers. At this point, 
it is useful to recall that the auxiliary vector A has only been introduced to get rid of the trace 
terms of the expectation values of the composite operators and has no physical significance. By 
undoing the contraction with A, these terms appear again. Consider as an example 

/KP,rn,2,l) = J0^...J^{A.q,nA.q,)f{ku...,k,p,m) (5.21) 

^ A^^A'^^A^^l^..J^q,,,,q,,,,q,,,J(k,,...,k^^^^^ (5.22) 

One notices that the way of distributing the indices in Eq. (5.22) is somewhat arbitrary, since due 
to the contraction with the totally symmetric tensor A'^'^ A^^ A^'^ , the result of the corresponding 
tensor integral can be taken to be fully symmetric as well. This is achieved by distributing the 
indices among the qi in all possible ways and dividing by the number of permutations one has 
used. Thus Eq. (5.22) becomes 

1 f d^ki f cPki 
Ii{p,m,2,l) = A'^^A'^^A'^^-y ^^...y ^^(?i,^,gi,^3?2,M, +gi,^,?i,^3?2,M2 

+(li,i^i(li,i^2(l2,,.3)f{ki,...,ki,p,m) . (5.23) 
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Generally speaking, the symmetrization of the tensor resulting from 



l[{A.q,r (5.24) 

1=1 

can be achieved by shuffling indices, [69-74] , and dividing by the number of terms. The shuffle 
product is given by 



(5.25) 



C {h, . .. , fei) LU (A;2, . .. , A;2)^ LU . . . LU {ki, . . . , fcj) 

m n2 ni 

where C is the normalization constant 

C=( ^ \\ (5.26) 
\ni,... ,nij 

As an example, the symmetrization of 

91,^2 52,/i3 (5.27) 

can be inferred from Eq. (5.23). After undoing the contraction with A in integral (5.18) and 
shuffling the indices, one may make the following Ansatz for the result, which follows from the 
necessity of complete symmetry in the Lorentz indices 

[N/2\+l j-1 N 
3=1 k=l l=2j-\ 

In the above equation, [ ] denotes the Gauss-bracket and {} symmetrization with respect to the 
indices enclosed and dividing by the number of terms, as outlined above. The first few terms 
are then given by 

i?o = 1 , (5.29) 
= ^iP/^i ) (5-30) 

^{miM2} = ^iPmiPm2 + ^zS-MiMa ■, (5-31) 
R{H1U.2HZ} = MPn^Pn^Pnz + ^2fl'{/.i/X2PM3} ■ (5-32) 

The scalars Aj have in general different mass dimensions. By contracting again with A, all trace 
terms vanish and one obtains 

Ii{p,m,n,,...,nj) = A^\..A'^^ R{^^,„^^y (5.33) 
= A.iA.pf (5.34) 

and thus the coefficient Ai in Eq. (5.28) gives the desired result. To obtain it, one constructs a 
different projector, which is made up only of the external momentum p and the metric tensor. 
By making a general Ansatz for this projector, applying it to Eq. (5.28) and demanding that the 
result shall be equal to Ai, the coefficients of the different Lorentz structures can be determined. 
The projector reads 

[iV/2]+l [N/2]-i+l N 

n....... = F{N) y: c{z,n)[ n ^^)( n (5-35) 

i=l 1=1 ^ k=2[N/2]-2i+3 ^ 
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For the overall pre-factors F{N) and the coefficients C{i,N), one has to distinguish between 
even and odd values of N, 

C^^(kN) = ( i)Nf2^^^U2 2""~''^'-'^"^(^ + + N/2 + k- 3/2) 

^ ' ^ ^ ^ r{N/2-k + 3/2)T{2k)T{D/2 + N/2-l/2) ' ^ ^ ^ 

^ ^^^^ - p - l)r(7V/2 + L>/2 - 1) ' ^^■'^^^ 
f;^7Vl - ( 1)^/^+^+^ 2^^-^/^-MA^ + l)r(g/2 + 7V/2-2 + fe) 

^ ^"'^^^ ~ ^ r{N/2-k + 2)r{2k-i)r{D/2 + N/2-1)' ^^■'^^^ 

F-^-(N) - 2^-^/^r(g/2 + 1/2) 

^ - (i? - l)r(7V/2 + D/2 - 1/2) ■ ^^■^^> 

The projector obeys the normalization condition 

U^„...,^^R^"-'^'' = A,, (5.40) 

which implies 

n,„„,,p^\..if-^i. (5.41) 

(5.42) 



As an example for the above procedure, wc consider the case = 3, 

.j^_3W^^p^2)^^^l^^^) . (5.43) 



^VlM2/i3 ~ D — 1 



Applying this term to (5.23) yields 

IiiP^rn,2,l) = (^^^e J J ^i-^P <i-<i^P-<i^ 

-p'^qlp.q2 + {D + 2){qi.p)^q2.p)f{ki,...,ki,p,m) . (5.44) 

It is important to keep p artificially off-shell until the end of the calculation. By construction, 
the overall result will not contain any term oc 1/p^, since the integral one starts with cannot 
contain such a term. Thus, at the end, these terms have to cancel, one can set p^ = and the 
remaining constant term in p^ is the desired result. 

The above projectors are similar to the harmonic projectors used in the MINCER-program, 
cf. [65,75]. These are, however, applied to the virtual forward Compton-amplitudc to determine 
the anomalous dimensions and the moments of the massless Wilson coefficients up to 3-loop 
order. 

The calculation was in general performed in Feynman gauge. For the external quark and 
gluon lines, the projectors (5.15, 5.13) are applied, which requires to include the ghost terms 
into the calculation. We also performed part of the calculation keeping the gauge parameter in 
i?^-gauges, in particular for the moments A^ = 2,4 in the singlet case and A^ = 1,2,3,4 in the 
non-singlet case, yielding agreement with the results being obtained using Feynman-gauge. In 
addition, for the moments A" = 2,4 in the terms with external gluons, we applied the physical 
projector in Eq. (5.14), which serves as another verification of our results. The computation of 
the more complicated diagrams was performed on various 32/64 Gb machines using FORM and 
for part of the calculation TFORM, [76], spending about 250 days of computational time. 
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5.3 Results 



We calculated the unrenormalized operator matrix elements treating the IPI-contributions ex- 
plicitly. They contribute to ^q^) ^55,(3 ^^'^ "^fl'Q^ ■ obtains the following representations 



+4>'(«).|^)n'"(o.|^). (5.45) 



Ca -c;-ft'"(o.f)-i::;n"'(o.^ 

-2iln«(o.f ) + i:>.(o,^)na.(o,f ) , (5.46) 

j(3),NS x(3),NS,irr . ^^2 

The self-energies are given in Eqs. (4.5, 4.6, 4.7, 4.9). The calculation of the one-particle ir- 
reducible 3-loop contributions is performed using MATAD. The amount of moments, which 
could be calculated, depended on the available computer resources w.r.t. memory and compu- 
tational time, as well as possible parallelization using TFORM. Increasing the Mellin moment 
by two demands both a factor of 6-8 larger memory and CPU time. We have calculated the 
even moments = 2, . . . , 10 for ^q^, ^^gg,Q-i ^'^^ ^^g!^' ^'-'^ ^Q?^^ up to = 12, and for 

.(3),NS 4(3),PS .(3) i Ar 1/1 

{i) Anomalous Dimensions : 

The pole terms of the unrenormalized OMEs emerging in the calculation agree with the general 
structure we presented in Eqs. (4.15, 4.24, 4.25, 4.37, 4.38, 4.45, 4.55). Using lower order 
renormalization coefficients and the constant terms of the 2-loop results, [29,33,35,86], allows to 
determine the fixed moments of the 2-loop anomalous dimensions and the contributions oc Tp 
of the 3-loop anomalous dimensions, cf. Appendix 8.2. All our results agree with the results of 
Refs. [18,19,23,24,79,80]. The anomalous dimensions 7gg'' and 7gg'*'^^ are obtained completely. 
The present calculation is fully independent both in the algorithms and codes compared to 
Refs. [18, 19, 23, 24, 80] and thus provides a stringent check on these results. 



{a) The constant terms a\- {N): 

The constant terms in Eq. (3.3) at O(a^), (4.15, 4.24, 4.25, 4.37, 4.38, 4.45, 4.55), are the new 
contributions to the non-logarithmic part of the 3-loop massive operator matrix elements, which 
can not be constructed by other renormalization constants calculated previously. They are given 
in Appendix 8.3. All other contributions to the heavy flavor Wilson coefficients in the region 
S> are known for general values of A^, cf. Sections 2, 4. The functions af^{N) still contain 
coefficients oc C2 and we will see below, under which circumstances these terms will contribute 
to the heavy flavor contributions to the deep-inelastic structure functions. The constant B4, 

^^Partial results of the calculation were presented in [77, 78]. 
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(4.10), emerges as in other massive single-scale calculations [81]. 

(m) Moments of the Constant Terms of the 3-loop Massive OMEs 

The logarithmic terms of the renormalized 3-loop massive OMEs are determined by known 
renormalization constants and can be inferred from Eqs. (4.17, 4.26, 4.27, 4.39, 4.40, 4.47, 
4.58). In the following, we consider as examples the non-logarithmic contributions to the second 
moments of the renormalized massive OMEs. We refer to coupling constant renormalization in 
the MS-scheme and compare the results performing the mass renormalization in the on-shell 
scheme (m) and the MS-scheme (m). 

For the matrix elements with external gluons, we obtain : 

{-'-^ + ?5?B. - 104a + - f C, + 64C. ln(2) 

^79 64 701 
+TfC'f I - y ^4 + 32C4 ^ ^Cs + 8OC2 - I28C2 ln(2) 

, 21586 3605^ \ 55672 889^ 128. 

+^^^^ I " W + 162"^^ j + ^^^^ + -81^' + T^' 

, 7054 704 \ f 22526 1024, 64, 



(5.48) 



,(3),MS. , __2 _ TC'f ^^^^^^ ^^B +72^ \+Tr r ( ^^^^^^ 

208^ 2330 \ ^^2/^ 8042 64^ 

+— B4 - 104C4 + -^Cs j + T,Cl i^-^ - -84 + 32C4 



81 / ' \^ 2187 ' 162 ' ^ " \^ 729 ' 81 

, 7054 704 \ A78 1024, , , , , 

+n,nc. n,nc, (- + — & I . (5.49) 

- „,r|C.|!^-^C3).n.r^C.(-^.i|ic3). (5.50) 

.<«(.^ = ™^2) = r.Ci(-^.fB.-72a.^C3) (5.51) 

^T.C.C. I ^ - ^B, + 104a - ?^C3 + - 64C. .n(2) 
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+TFCj^ ( W + " ^^^^ + "^^^ " + ^^^^^ ^''^^^) 



^^2^ 21586 3605 \ 55672 889 128 

" l62"^^J + ^^^^ |^^29~ " " T'^' 

/ 57226 1408 \ / 29908 2048, 64, 



(5.52) 



,(3),Ms. o _ o 2X _ ^ ^2 f 174055 88 29431 



rr. ^ ^ /123II3 208^ _ 2330 \ _ . ^8042 64^ 



81 \^ 2187 162 ^ \^ 729 81 

57226 1408^ \ /^6904 2048 \ , 

Comparing the operator matrix elements in case of the on-shell scheme and MS-scheme, 
one notices that the terms ln(2)(^2, C2 ^-re absent in the latter. The C2 terms, which contribute 
to a\^\N), are canceled by other contributions through renormalization. Although the present 
process is massive, this observation resembles the known result that (^2^terms do not contribute 
in space-like massless higher order calculations in even dimensions, [82]. This behaviour is 
found for all calculated moments. In addition, (^4-terms occur, which may partly cancel with 
those in the 3-loop hght Wilson coefficients, [26]. Note, that Eq. (5.50) is not sensitive to mass 
renormalization due to the structure of the contributing diagrams. 

An additional check is provided by the sum rule, [31], 

A'SliN ^ 2) + Af^^QiN = 2) + A^^^q{N = 2) = 0, (5.54) 

which is fulfilled in all renormalization schemes and as well as on the unrenormalized level. 

Unlike the operator matrix element with external gluons, the second moments of the 
quarkonic OMEs emerge for the first time at O(a^). To 3-loop order, the quarkonic OMEs 
do not contain terms oc ^2- Due to their simpler structure, mass renormalization in the on-shell- 
scheme does not give rise to terms cx C2, ln(2)C(2). Only the rational contribution in the color 
factor oc TpCp turns out to be different and A^^'q \ (5.56), is not affected at all. This holds 
again for all moments we calculated. The non-logarithmic contributions are given by 

,(3),MS,PS. 2 _ ^2 ^ ^ 830 64 1280 V ^ ^2 95638 

^Q? (/^ -"^,2j - ^P^^^^l ^^87 "9" ^ ~ ^^^^^^ J ^29~ 

9 ^ 81 ^ / ^ I 2187 81 ^ 
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, 34312 1024^ 



81 ' \^ 2187 81 

= -^ 2) = n,nc, (-^ + . (5.56) 



128 „ 848. \ / 25024 1792 



B4 + 64C4-^C3 +r/C'FhTT^-^C3 



9 ' 9 ' / ^ \ 2187 81 

-ipB,.64a-^.).r^cif^-lZfc. 



/ 46336 1024 



256^ 17168, \ / 26056 1792, 

+-B,-128a+— Csj + — Cs 

^„,T|cil|I?-15|!c3|. (5.59) 



729 27 



.(3)MS. 2 X / 101114 128^ , 8296 \ 436094 

2) = r.cpc. _ - —B, + 128a - ^c, + t,ci -— 
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Finally, the sum rule, [31], 

AgfS(Ar = 2)+A2^^^(Ar = 2)+>lX(Ar = 2)+>l2Q(7V = 2) = (5.61) 

holds on the unrenormalized level, as well as for the renormalized expressions in all schemes 
considered. 

FORM-codes for the constant terms a^y (iV), Appendix 8.3, and the corresponding moments 
of the renormalized massive operator matrix elements, both for the mass renormalization carried 
out in the on-shell- and MS-scheme, are attached to this paper and can be obtained upon 
request. Phenomenological studies of the 3-loop heavy flavor Wilson coefficients in the region 
Q'^ » will be given elsewhere [83] . 



6 Heavy Quark Parton Densities 

In the kinematic region in which the factorization relation (2.5) holds, one may redefine the 
results obtained in the fixed flavor number scheme, which allows for a partonic description at 
the level of (n/ + 1) flavors. As before, we consider Uf massless and one heavy quark flavor. 
Since parton distributions are process independent quantities, we define the parton distributions 
for (n/ + 1) flavors from the light-flavor parton distribution functions and the massive opera- 
tor matrix elements for Uf light flavors. Also in case of the structure functions associated to 
transverse virtual gauge boson polarizations, like F2{x,Q'^), the factorization (2.5) only occurs 
far above threshold, ~ Arn^x/ (1 — x), and at even larger scales for Fl{x, Q^). The following 
set of parton densities is obtained, cf. [31] : 

hirif + 1, /x^ m^ N) + fj^irif + 1, //^ m^ N) = (^rif, ^, ■ [Muf, N) 

+^^Iq {rif,^,Nyj:{nf,t,',N) 

+A,,,Q (nf,-^,NyG{nf,pi^N), (6.1) 

/Q(r^/ + l,//^m^A^) + /Q(n^ + l,/x^m^A^) = Aj^ (^n^, ^, AT^ • E(n/, /x^, AT) 

+AQ,(^nf,^,N^ .G{nf,i,',N) . (6.2) 

Here, fk{fk) denote the light quark and anti-quark densities, fqifQ) the heavy quark densities, 
and G is the gluon density. The flavor singlet, non-singlet and gluon densities for (n/ + l) flavors 
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are given by 



A 



NS 

qqA 



+^Qq 



,2 



+ 



m 



^f^qg^Q r^/' — ^' ^ + ^Q9 — ^ 



A{nf + l,f/,m^,N) = 



fkiuf + 1, //^ TV) + f^{nf + 1, m^, N) 
^ -E(n/ + l,/x^m^7V) 



nf + 



(6.3) 
(6.4) 



(6.5) 



Note, that the new parton densities depend on the renormahzed heavy quark mass m^. As 
outhned above, the corresponding relations for the operator matrix elements depend on the 
mass-renormahzation scheme. Furthermore, m — m{as{n'^)). This has to be taken into account 
in QCD-analyzes, in particular cannot be chosen constant. 

The normalization of the quarkonic and gluonic operators obtained in the light-cone expan- 
sion can be chosen arbitrarily. It is, however, convenient to select the relative factor such, that 
the non-perturbative nucleon-state expectation values, S(n/,/x^,iV) and G(n/, /x^, iV), obey 



E(n/, A^^ = 2) + ^(n/, A^^ = 2) = 1 



(6.6) 



due to 4-momentum conservation. As a consequence, the OMEs fulfill the relations (5.54, 5.61). 
The parton densities (6.1-6.5) can be applied in other hard-scattering reactions at high energy 
colliders in kinematic regions where the corresponding power corrections oc (m^/Q^)'^, A; > 1 
can also be safely disregarded. 

Conversely, one may extend the kinematic regime for deep-inelastic scattering to define the 
distribution functions (6.1-6.5) upon knowing the power corrections which occur in the heavy 
flavor Wilson coefficients Hl{x,Q'^ / (j? ,111^ / ^^). This is the case for 2-loop order. We separate 



H 



j.asymp 



Q ^ \ I ii:?, power 
9"' 9" ) i 



9 ' 9 



9 ' 9 



(6.7) 



where Hf^^^'"''(x, Q^Z/i^, m^Z/U^) denotes the part of the Wilson coefficient given in Eq. (2.5). If 
one accounts for H\''^"^^\x,Q'^ / [j? / ix^) in the fixed flavor number scheme, Eqs. (6.1-6.5) are 
still valid, but they do not necessarily yield the dominant contributions. In the region closer to 
threshold, the kinematics of heavy quarks is by far not collinear, which is the main reason that 
a partonic description has to fail. Moreover, relation Eq. (1.4) may be violated. In any case, 
it is not possible to use the partonic description (6.1-6.5) alone for other hard processes in a 
kinematic domain with significant power corrections. 
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7 Conclusions 



We calculated the 3-loop massive operator matrix elements, which form the heavy flavor Wil- 
son coefficients, (2.11-2.15), together with the known massless Wilson coefficients in the region 

3> rv? due to the factorization theorem (2.5). All but the power-suppressed contributions are 
obtained in this way. Furthermore, all operator matrix elements needed to derive massive quark- 
distributions at the 3-loop level were calculated. We presented in detail the renormalization of 
the massive operator matrix elements, leading to an intermediary representation in a defined 
MOM-scheme. This is necessary to maintain the partonic description required for the factoriza- 
tion of the heavy flavor Wilson coefficients into OMEs and the light flavor Wilson coefficients. 
The representation of the heavy flavor Wilson coefficients in the asymptotic region, effectively 
reached for the structure function ^2(^,(5^) for Q'^jvr? ~ 10, is available for general values of 
N in analytic form, up to the constant parts af^{N) of the unrenormalized 3-loop OMEs. A 
number of fixed values of Mellin moments for these constant parts were calculated, reaching 
up to A^ = 10, 12, 14, depending on the complexity of the corresponding operator matrix element. 
Although general methods are available to reconstruct the recurrence formulae for anomalous 
dimensions and Wilson-coefficients as a function of A" by a finite number of moments, [84], the 
number of moments calculated for a\f{N) is still far too low. Through the renormalization of 
the massive OMEs, the corresponding moments of the complete 2-loop anomalous dimensions 
and the Tp-terms of the 3-loop anomalous dimensions are obtained, as are the moments of the 
complete anomalous dimensions 7gg'*'^^(A^) and 'yql\N), which agree with the literature. 

The results were presented performing the coupling constant renormalization of the OMEs 
in the MS-scheme and the mass renormalization in the on-shell scheme. After a transformation 
to the MS-mass, the (^2^terms are canceled completely. Although being a massive calculation, 
which is indicated by the emergence of the number B4, the use of the MS-scheme moves the 
structure of the result towards those observed in massless 3-loop calculations. 
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8 Appendix 



8.1 Feynman Rules 

For the Feynman rules of QCD, we follow the convention of Ref. [55]. D-dimensional momenta 
are denoted by pi and Lorentz indices by Greek letters. Color indices are denoted by a, 6, and 
z,j are indices of the color matrices. Solid lines represent fermions and curly lines gluons. The 
Feynman rules for the quarkonic composite operators are given in Figure 8. Up to 0{g^) they 



N-l 



N > 1 



P2,J 





\N-l-2 



g^A'^A^j^ Epo' T.l-U^P2V{Ap, 
TV > 3 




P5,P,c 



ff3A,A.A,A7± EZzU^.poyi^.Pir-"--' 

[ {tH''t'')ji{A.p4 + A.ps + A.pi)'-^-i(A.p5 + A.pi)™-'-i 



{m'')ji{A.pi + A.p5 - 


h A.pi) 


'-i-^A.p4 + A.pi)"^- 


-1 


{tHH^)^i{A.p-i + A.p5 - 


h A.pi) 


'-J-i(A.p5 + A.pi)™- 


-1 


{tHH'')ji{A.p3 + A.p5 - 


h A.pi) 


'-J-i(A.p3 + A.pir- 


-1 


{tHH'>)ji{A.p3 + A.pi - 


h A.pi) 


'-^-\A.pA + A.pi)"'- 


-1 


{tHH'')ji{A.p3 + A.pi - 


h A.pi) 


'-J-i(A.p3 + A.pir- 


-1 



7+ = 1 



iV > 4 
75 • 



Figure 8: Feynman rules for quarkonic composite operators. A denotes a light-like 4-vector, 
is an integer. 



can be found in Refs. [10] and [85]. Note that the 0{g) term in the former reference contains a 
typographical error. We have checked these terms and agree up to normalization factors, which 
may be partly due to a different convention in the standard Feynman rules. We newly derived 
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the rule with three external gluons. The terms 7± refer to the unpolarized (+) and polarized 
(— ) calculation, respectively. Gluon momenta are taken to be incoming. The Feynman rules 
for the unpolarized gluonic composite operators are given in Figure 9. Up to 0{g'^), they can 
be found in Refs. [11] and [16]. We have checked these terms and agree up to 0{g^). At 0{g), 
we agree with [11], but not with [16] and [55]. At 0{g^), we do not agree with either of these 
results, which even differ from each other 



p, V, b p, a 



g^,(A • pf - (A^p, + A,p^)A ■ p + p'^A^A, 



N >2 




{A,gx^ - Ax9^^)A ■ pi + A^(pi,,Aa - Pi,xA,)\ (A • pif-^ 
A ■ piPa^pA^ + A • p2Pi,t,A^ - A ■ piA ■ p^g^^v - Pi ■ P2A^A,, 



Ef=o'(-A-PiV(A-p2 



j P1^P2^P3^P1 l I j P1^P3^P2^P1 \ 




g2i±ur \^r''^r'i^o,^,„ipuP,,p,,P4) 

+ r'.f'''0,x.APl,P3,P2,P.i) + /""^/''^"OM-A(Pl, P4, P2, P3) 
O^^XaiPl,P2,P'i,P4) = A^Aa< -5^^(A -ps + A-p4,)^-'^ 



-[P4,mA. - A • p45m.] ES'(A • P3 + a • p,yiA . p,r-^-^ 



-bi,.A, - A-p,g,,] j:l-\-A.p,y{A-p, + A-p,r-^-^ 
-[A ■ piA ■ pig^„ + pi ■ PiA^A„ - A ■ piPi^„A^, - A ■ piPi^^,A„] 
X Ef=o' E ■=o(-A • Pif-'-'iA • P3 + A . p4)^--'(A • P4>'j 

1 jl^U j 1 A*— !-(T J 1 fl-^f, X<-^<T j ' — 



Figure 9: Feynman rules for gluonic composite operators. A denotes a light-like 4-vector, = 0; 
is an integer. 



^^We would like to thank J. Smith for the possibility to compare with their FORM-code used in Refs. [29,56, 
86,87], to which we agree. 
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8.2 The 3— loop Anomalous Dimensions 

The 3-loop anomalous dimensions Jqq{N) and JqgiN) and the contributions oc Tp to 
Iqq^^ (N) , 'ygq{N) and '~igg{N) are obtained from the single pole terms in the present calcula- 
tion for even values of N and for \^^^{N) for odd values of N. In the latter case, also ^q'^^{N) 
with 7q^'^^(-/V) = 'jq^^^iN) +^q'^^{N) cShU be obtained, which will be considered elsewhere [83]. 
The anomalous dimensions are : 



TpCp 



^ -^^ 243 3 V ^ ^/^•^ 243 ' 



14728 



^ ,^ 618673 968 /^ ^ \^ 2485097^ 
+ '"^)^-151875 + 7^ (^^ - + ^^^^ 



506250 



(8.1) 



2217031 



675000 



, ^ ,^ 126223052 3872 
-(l + 2n/)T^ 72930375 + 



1988624681 
^ ^ 4084101000 ^ 



(8.2) 



11602048711 
^10210252500 ^ 



F^F 



,^ ^ ,^ 13131081443 2738 / ^ \^ 

1 + 2nt)Tp \ [Cf-Ca]^ 

' ^' -^13502538000 945 V 



343218329803 
648121824000 ■ 



(8.3) 



39929737384469 
^22684263840000 ^ 



/, ^ 265847305072 50176 /^ ^ \ ^ 
^ ^' 420260754375 27225V ^ ^)^' 



A 1,3 



(8.4) 



1028766412107043 839864254987192 



1294403123475000 



485401171303125 



2566080055386457 49928 / 

-(1 + 2n^)TF \ h 

^ ^' 5703275664286200 39039^ 



[pF — C'a^ Cs 



(8.5) 



69697489543846494691 86033255402443256197^ 



83039693672007072000 



54806197823524667520 



^.6) 
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7j?(4) 



8464^ 1384 



243 
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Cf) + y (^-416CaC^ + 288C2 + 128C2 ^ 



7178 2 ,556 8620 2 



(8.7) 



/4481539^ 9613841 ^ 
+ ^^^^^H ^750"^^ + 3037500^^; 



^ ( 2832C;^ - 3876CaCf 



+1044^1) 
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278546497 
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H — CaC 
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(1 + 2n;)7V (yY^^^a ^ 340341750 
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24648658224523 
289340100000 

4374484944665803 
56710659600000 
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10125 15 
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agree with the anomalous dimensions given in [18, 19,22-24]. 
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^(3) 

8.3 The O(£0) contributions to A^. (AT) 

The constant contributions to the unrenormahzed massive operator matrix elements at 0{a^) 
read : 
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2903694979 \ 2 4936013830140976263563 _ 25088 12544 

^ 6670805625 + f f 80731922811135750000 ~ 81675 ^ 9075 

94499430133 282148432 \ 2^ 430570223624411 802816 
^ 1940598000 ^ 4002483375^^ J + f f 2780024890190625 ~ 735075^^ 

319072 \ . f 454721266324013 175616 

^11026125^^ ] + ^/ F "624087220246875 ^ 735075^^ 

547424 ^ , , ^ , 

-24257475^' ' '^'^^^ 

T C C ( 968307050156826905398206547 12482 ^ 12482 
^ ^ I 107727062441920086477312000 ^ 117117 ^ ~ 13013 



64839185833913, 489403711559293 , 

^C3 + -.oon.-.nnonnr-.nnn C2 



16206444334080" 1382612282251200 

T (72 I _ 190211298439834685159055148289 _ 24964 ^ 12482 
^ ^ ' 2962494217152802378126080000 ~ 117117 ^ 13013^^ 

418408135384633 72904483229177 
^ 8103222167040 ~ 15208735104763200 

2^ ( 1727596215111011341 _ 798848 11471393 
^ ^ ^ \^ 13550982978344011200 ~ 1054053^^ ^ 347837490^^ 

f 6621557709293056160177 24964 ^ 
^ I 12331394510293050192000 150579 

^''"''"^^C. I . (8.42) 



63306423180 



(3),ps... / 100096 896 256 . 

<Q (2) = rifTpCp I + -^Cs - -^C2 ) (8.43) 



(3)PS,,x / 118992563 3388^ 4739 ^ , 
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(3),PS/^x 



%q,Q yp) 



aSf (10) = n^ncp 



2/ 



«Sf (14) = n,T|C^ 



17732294117 1936 _ 9794 
"10210252500 ^ 2835^^ ~ 694575^^^ 

20110404913057 1369 135077 
"27221116608000 ^ 3645 ^ 4762800 y 

308802524517334 175616 4492016 
"873722108345625 ^ 735075 ^ 121287375 

6724380801633998071 24964 583767694 
"38535607844665781850 ^ 150579^^ ^ 15826605795^^ 

616164615443256347333 22472 189601441 
" 7545433703850642600000 ^ 184275^^ ^ 5533778250 y 



(8.45) 
(8.46) 
(8.47) 

(8.48) 



(8.49) 



(iii) ag](iV) 



«S(2) 



_ , 170227 88^ „ , 31367, 1076, \ ^ ^ ^ / 154643 



208^ 7166, \ ^ ^2 / 15574 64^ 

+— B4 - 104C4 + - 54C2 + TfC], - -B4 + 32C4 



3421 704 \ rT.2r^ 20542 4837 670 



11696 
729 



569 . 256 



64, 



6706 616 



250 



9 



27 



2187 81 81 



, 158 896^ 40^ 



(8.50) 



«S(4) 



2 ( 425013969083 559 2124 352717109 4403923 
^ ^ ' 2916000000 ~ ^ ^ ^ ^5~^^ ~ 5184000 ~ 270000 ' 



^ ^ ^ , 95898493099 646 2907 172472027 923197 

+TfCfCa I cr,An^r^r^,^r, + "?T^^4 ^^C4 H o^^r.r.r. Ca 7^T77:7rC2 



874800000 25 



25 



864000 



40500 



_ ^2 / 87901205453 174^ 783, 937829, 62019319, 



699840000 25 



25 



12960 



3240000 



( 960227179 1873781 . 120721 



+ 



\ 29160000 51840 



Ca + 



13500 



■C2 I +T^Cf\ 



1337115617 
874800000 



73861 , 8879111 



324000 



810000 



176 3 , ( 947836283 18172 
:-/ i^Ca + nfl pUA I ...^ Ca 



135 



72900000 2025 
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11369 \ 2^ / 8164734347 130207 1694939 , 



13500^ / ■ J ^ ^ U374000000 20250 810000 

T /_48989733311629681 _ 2938^ 17466 _ 748603616077 
^ ^ I 263473523712000 ~ "315" ^ 245 ~ 11379916800 

_ 93013721 \ / 712876107019 47332 23666 

~ 3457440 J ^ ^ ^ ^ I 55319040000 ^ 2205 ~ 245 

276158927731 4846249 \ ^ ^2 3873986781136411; 
^ 1896652800 ^ 11113200 ^^J ^ 137225793600000 

2480 1240 148514798653 4298936309 ' 
~ 441 ^ "^ ^9~^^ ^ 711244800 ^ 388962000 

/ 706058069789557 3393002903, 6117389, 
^ ^ I 18819537408000 116121600 555660 

/ 447496496568703 666922481 49571129, 
^ ^\ 54890317440000 284497920^^ 9724050 ^^ 

176 3 2^ f 12648331693 4433 23311 ^ 

" 189^^^' + ""^^^^-^ 1^ 735138180 " '567^' + im32^^^ 

f 8963002169173 111848, 11873563, A 

^nfTlCF + Cs + C2 (8.52 

y 1715322420000 19845 19448100 / ^ ' 

rp (j2 ( 358497428780844484961 899327^ 64021 
^ ^ I 2389236291993600000 ~ 113400 ^ 1050 

_ 12321174818444641 _ 19581298057 \ ^ ^ ^941315502886297276939 
112368549888000 ~ 612360000 + f f a y §352327021977600000 

515201 g _ 515201 5580970944338269 495290785657 ' 
^ 28350 ^ ~ 6300 ^ 56184274944000 ^ 34292160000 

T ^2 23928053971795796451443 _ 749^ 749 719875828314061 
^ ^ \^ 36585180721152000000 162 ^ "36"^^ "*" 1404606873600 

2484799653079 , \ ^2^^56313300657148129 58802880439, 
480090240000 j ^ ^ I 4147979673600000 2388787200 

46224083 \ ^.2^ ( 986505627362913047 185046016777 
^ 4082400 / ^ ^ ^ I ~ 87107573145600000 ~ 50164531200 



7527074663 \ 296 3 ^ f 24718362393463 125356 

3429216000^^ ~ 405 ^^^^'^^ ^ ^ 1 1322697600000 ~ 18225 
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2118187 \ 2^ / 291376419801571603 887741 

^2916000^^ I ^^-^ ^ ' ~ 32665339929600000 ^ 174960^^ 

139731073 , , 

C2 (8.53) 



1143072000 

(3) _ 2/^6830363463566924692253659 563692 483988 
^^^^ ' ~ ^ ^ y 685850575063965696000000 ~ 81675 9075 

103652031822049723 20114890664357 

"Ca roi -.mnrvAAAA C2 



415451499724800 581101290000 

'872201479486471797889957487 1286792 643396 
^ ^ P ^\ 2992802509370032128000000 ^ 81675 ^ ~ 9075 

761897167477437907 , 15455008277 , 

-C3 + ...o.^o^r C2 



33236119977984000 660342375 

_ ^2 / 247930147349635960148869654541 11808^ 53136, 
^ ^ ' 148143724213816590336000000 3025 ^ 3025 

9636017147214304991 14699237127551 
^ 7122025709568000 ^ 15689734830000^^ 

,'23231189758106199645229 123553074914173, 4206955789, 
^ ^ ' 633397356480430080000 5755172290560 377338500 

^2^ / _ 18319931182630444611912149 _ 502987059528463 
^ ^ ^ ' ~ 1410892611560158003200000 ~ 113048027136000^^ 

24683221051 \ 896 ^3 ( 297277185134077151 

^ 46695639375 ~ 1485 + ^•'^ ^ ^ 15532837481700000 

1505896 189965849 A ^ ( 1178560772273339822317 
<3 + iQQ««no^n^2 + n/ij^Cp 



245025 ^ 188669250^ y ^ - ^ 107642563748181000000 

62292104 49652772817 . 
^13476375^^ ~ 93391278750^^ ' ' ^ ' ' 



(3) 



(3) ^2^ /83204 616, 290, \ / 5000 896, 4, 



.55) 
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(3) _ / 835586311 18172 71899 

21270478523 130207 _ 1401259 
^ ^ ' 874800000 ^ 20250 ~ 810000 



(3) ta^ _ / 277835781053 4433, , 2368823, 

%,q(6) - njT^CA\ 5881105440 " + ^55660"^^ 



2^ I 36123762156197 111848 26095211 
^^•'^ ^ ^ ' ~ 1715322420000 ^ 19845 ~ 19448100^^ 

(3) C8^ - T'^c ( 15''327027056457 125356 7917377 
%g,Q\ ) - F a\ 3958092800000 18225 2268000^^ 

^2^ / 201046808090490443 887741 3712611349 
'^^f ^ ^1 ~ 10888446643200000 ^ 174960^^ ~ 3429216000^' 

(3) _ ^2^ ( 6542127929072987 _ 1505896 1109186999 

%9,Qy ) - F A y 191753425700000 ~ 245025 ^ 377338500 

^.^^ / 353813854966442889041 62292104, 

+nfTiCF\ \ Cs 

' 21528512749636200000 13476375 

83961181063 

;C2 



93391278750 



(3) ^21 - rrr f ^^^^^^ ^^^b +i28r ^^^^r 



-^^^^ (- w - f - - If ^= 1 + ^^^^ I 

2048, 416, \ / 92200 896, 208, 

+ ^^^J + ^^^^^^ [-729- - Tf^' + 

J3) .4^ _ 7. 1^ 5501493631 176 176 8258 13229 

S..q(4) - TfCfCa ^ - 218700000 " 4^^^ + " lo^^^ + "SlOO"^^ 

^^2/ 12907539571 352^ 176, 132232. 398243^ 

+iFCp 1 Da C4 H Cs C2 

^ ' 145800000 45 5 2025 27000 

2^ 1914197 2816 1252 
(^"911250" + l3^^' + "67^^^ 
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/ 50305997 1232, 626, , 

^ ^ ^ , 384762916141 704 704 240092 403931 

-/pCpG^ I Da -\ Ca — Ca H C2 

' 24504606000 315 35 19845 463050 

2 ( _ 40601579774533 1408 _ 704 27512264 
^ ^ ' 918922725000 ^ "SlS" ^ ~ "35"^^ ^ 694575 

24558841 \ 279734446 11264 8816 

' 3472875 ^^J ^ ^ 364651875 ^ 945 ^^^ 33075^^ 

,'4894696577 704, 4408 , , 

. 10318865954633473 _ 296 g 296 1561762 30677543 
^ ^ ^ ' 816633498240000 ~ 189 ^ ^ '21^^ ~ 178605 ^ 85730400^^ 

^^2 1 305405135103422947 592^ 296, 124296743, 

+TfCI + B4 C4 + Cs 

11432868975360000 189 21 ^ 4286520 

4826251837 .A „ f 864658160833 4736 . 12613 



1200225600" / ^ V 567106596000 567 ' 59535 

^2^/^9330164983967 296, 12613 . \ ,^ 
+''f^^^^ 1^1134213192000 " "sT^^ " 119070^^ J ^^"^^^ 

T C C \ 1453920909405842897 1792 ^ 1792 1016096 
^ ^ ^ ' 130475834846280000 ~ 1485 ^ ^ 165"^^ ~ 147015 

871711 \ ^ 2 11703382372448370173 3584 ^ _ 1792 
^ 26952750 + f f 667205973645750000 ^ 1485 ^ ~ 165'''^ 

62282416 6202346032 \ ( 1346754066466 28672 

^ 2695275 ~ 2547034875 + f f 756459357875 + 4455 

297472 \ 2^ / 4251 185859247 12544 148736 , 

"735075^^ J + ""^^^^^ [ 756469357875 " l45^^^ ~ 735075^^ ' ^^'^^^ 

T C C \ 1515875996003171876913331 _ 1264^ 1264 _ 999900989 
^ ^ ^ ' ~ 147976734123516602304000 ~ 1287 143 173918745^^ 

693594486209 . \ m r^2 I 48679935129017185612582919 2528 _ 



3798385390800' y ^ \ 4069360188396706563360000 1287 

1264, 43693776149, 2486481253717, 



143 ' 2260943685 ' 1671289571952 
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2^ I 2105210836073143063 20224 _ 28514494 
^ ^ ' "1129248581528667600 ^ 3861 ~ 57972915^^ 

ri T'^-C (^9228836319135394697 _ 8848 14257247 . 

^' ^ I 2258497163057335200 ~ 386l''^ ^ 57972915^^ ' ^ ' ' 

(3) - T C C ( 1918253569538142572718209 3392 3392 

%g,Q[ > - F F A I - 199199449731556964640000 ~ 4095 * ^ l55"^' 

2735193382 1689839813797 

"Cs - riior,iiiAOAAAC2 



553377825 5113211103000 

2 / 143797180510035170802620917 6784 3392 
^ ^ ' ~ 17429951855894984406000000 ^ 4095 ^ ~ 455 

12917466836 4139063104013 
^ 774728955 ~ 4747981738500^^ 

, 337392441268078561 54272^ 98112488 , 
^ ^' 179653183425015300 12285^^ 184459275^^ 

2^ / 222188365726202803 3392 49056244 
^^•'^ ^ ^ 1 71861273370006120 ~ 1755^^ ~ 184459275^^ ' ' ^ ' ' 

f3) / X ^ ^9/^ 170227 88^ „ , 31367^ 1076 \ ^ ^ ^ / 154643 



208„ 7166 \ „ ^2 / 15574 64. 



B4 + 104C4 - -^(3 + 54C2 + TpC^ — — + — B4 - 32C4 



9 ' 27 ' / \ 243 9 

3421 704 \,^2r ( 20542 4837 670 V ( 11696 



569 256 \ 64 3 ^ I 76498 1232 40 

-n.T|C.|f-lZfc3-fG) (8.67) 

(3) _ ^2 f 29043652079 533 4698 610035727 92341 
S.,q(4) - T^C^y 291600000 2592000 + "6750"^' 

T C C ( 272542528639 1088 ^ 4896 3642403 73274237 
^ ^ ^ ^ 1 874800000 25~ ^ ^5~^^ ~ 17280 ^ 810000 
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2 ( 41753961371 44 _ 198 2676077 _ 4587077 
^ ^ ^ 1 1749600000 ^ 25 ^ ~ '25^^ ^ 64800 ~ 1620000^^ 

T^C I _ 1192238291 _ 2134741 _ 16091 \ /^_785934527 
^' 14580000 25920 ~ 675 J ^ ^ ^ I ~ 43740000 

32071 226583 \ 64 . ^ ( 271955197 13216 



8100' 8100 ' / 27 J ^ -1 1822500 405 

6526, \ , / 465904519 6776, 61352 \ 

""67^^'J + ""^^^^^ [~ 27337500 " 2025^^ " 10125^'J 

^^2/^37541473421359 56816^ 56376, 926445489353, 11108521 
^ 448084224000 2205 ^ 245 ^ 2844979200 ^ 555660 

/ 18181142251969309 _ 114512 57256 _ 12335744909 
^ ^ ^ ^1 54890317440000 2205 245 67737600 

94031857 \ ^ 2/^16053159907363 352 ^ 176 3378458681 
^ 864360 ^^j ^ ^' ^' I 635304600000 ^ 441 ^ Ig"^^ ^ 88905600 

8325229 \ 2^ /^ 670098465769 25725061 96697 
~10804500^y ^ ^ ^1 6001128000 259200 ~ 2835 

^2^ / 8892517283287 12688649 2205188 
^ ^ ^ \~ 490092120000 ~ 2540160 ~ 77175 

64 3 , [ 245918019913 3224 250094 

+^^fC3 + Ufl^CA y- 1312746750 + ~ 19845^' 

/_ 71886272797 3872 496022 \ 
+n/ ^ ^ 3403417500 ~ 2835^^ ~ 77175 ^ ' 

_ ^2 /^512903304712347607 108823^ 162587, 2735007975361, 
^ 18665908531200000 3780 ^ 630 6502809600 

180224911, \ ^ ^ ^ /^ 13489584043443319991 163882^ 81941. 

H C2 +-/fCfG4I D4 H Ca 

7654500 ^ y ^ ^ 43553786572800000 2835 315 ^ 

_3504113623243 414844703639 \ j: ^2 ( 5990127272073225467 
25082265600 ^ 3429216000 J + f I 228657379507200000 

37 37 3222019505879 12144008761 
^81 ~ 18^^ ^ 87787929600 ~ 48009024000^^ 

2^ I 16278325750483243 871607413 591287 
^ ^ ^ ' 124439390208000 7962624 ~ 14580 
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2^ I 7458367007740639 _ 291343229 _ 2473768763 
^ ^ ^ ' ~ 408316749120000 ~ 52254720 ~ 85730400 

64^3 /_ 102747532985051 54208 737087 

^27 + ^/ ^ I 486091368000 ^ 1215 ~ 51030 

2^ ( 1145917332616927 2738 70128089 . 
^ ^1 51039593640000 ~ 3645^^ ~ 10716300^^ ' ^ ' ' 



(3) / ^^9/ 15434483462331661005275759 17788828^ 17746492 

Q flOl = T C I I B ( 

^^'^^ ' 327337774462347264000000 571725 ^ 63525 

26909447654952 1109 1444408720649 
^ 519314374656000 ^ 55468759500 

T C C i 207095356146239371087405921 _ 35662328^ 17831164 
^ ^ ^ 771581896946961408000000 571725 ^ ^ 63525 



3288460968359099 , 6078270984602 ^ 



37093883904000 ' 46695639375 



^ ^2 / 553777925867720521493231 896 ^ 4032 , 

+rFCp ^ B4 C4 

20667372239650752000000 3025 3025^ 



7140954579599 _ 282148432 
' 198717235200 ~ 4002483375 



63059843481895502807 85188238297 . 33330316 



433789788579840000 729907200 735075 



T'^C ( _ 655690580559958774157 71350574183 _ 3517889264 
^ ^ ^ ' ~ 35787657557836800000 ~ 12043468800^^ ~ 121287375 

64^3 ^2^ 6069333056458984 215128 81362132 

^27 ^^^^^f ^ 26476427525625 ^ 4455 ~ 5145525 

^2^ / 100698363899844296 351232 799867252 . 
^ ^ ' 4368610541728125 ~ 735075^^ ~ 121287375^^ ' ' ^ ' ^ 



, ... (3),NS/,.N 



= (8-72) 

%-(2) = T.C.cJ^ + ^B.-64a+l^C3-^C.V^-^^f^ 
,Q \ ) ^ ""1 2187 9 ^ 81 81 / -^l 2187 
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- — B4 + 64(4 - -g-Ca + -gpC2j + T^C^F (^-^ - -gpCa 
512. \ ^2^/^ 100096 896 256 , 

■^C2 j + %r^c^ 1^-^^ + — C3 - ^C2 1 (8.73) 

I + l^B. - lOOC. + l^Cs + 

' 34992 9 ^162 ^^162^ 

_ ^0 , 35091701 200^ 1315, 29035, 

I ^3996^ - + '''^^ - -^^^ + 

, ( 188747 3200 830 
~^748 ST^^"^^' 

, 1271507 1400^ 415^ , 

1 -^749^ + - ) ^^-^^^ 

^ ^ ^ , 419369407 628^ 628 , 515597 . 10703 , 

JfGfG/i I 1 D4 — Ca H C'^ H C2 

^ ^ ^ ' 21870000 45 5 ^ 4050 ^ 4050 ^ 

2/^137067007129 1256 628 41131 4526303 
1^ 437400000 4^^^ + ~ + ^lOOO"^' 

T'^C I _ 1^1928299 _ 20096 _ 26542 
^ ^ ^ ' 5467500 405~^^ ~ 2025 

2^ ( 1006358899 8792 13271 , 

I - 10935000 + lo^^^ - -2025 • ^^-''^ 

^ ^ ^ , 816716669 728 _ 728 , 12569 , 16103 , 

I pCjpCj A I 1 D4 — Ca H Cfi H C2 

' 43740000 45 5 ^ 81 ^ 4050 ^ 

2/13213297537 1456 728 142678 48391 
1^ 36450000 4^^^ + ~ ~&J^^' + ^50~^^ 

I 9943403 23296 31132 
^ ^ ' ~ 303750 ~ 405 " 2025 

2^ / 195474809 10192 15566 

^ ^ ^ , 1541550898907 5672^ 5672 . 720065 , 1016543 , 

JfGfG/i I 1 D4 — Ca H C'^ H C2 

' 105019740000 315 35 3969 198450 

2 186569400917 _ 11344 5672 _ 7766854 55284811 
^ ^ ^ I 463050000 315~ ^ ^5~^^ ~ 33075 ^ 771750 
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2^ J 26884517771 _ 181504 _ 1712476 
^ ^ ' 729303750 2835 ~ 99225 

2^, 524427335513 11344 856238 

I ~ 4375822500 + " "9922^^^ ' ^^"^^^ 

/ 5307760084631 ^^^g _ ^1^^ 781237 19460531 
^ ^ ^ I 672126336000 ^ los" ~ ~35~^^ ^ 3780 ^ 3175200 

2/^4900454072126579 4108 6162 8425379 
^ 11202105600000 lOS" ^ ^ ~35~^^ ~ 33075 

1918429937 A . f 8488157192423 65728 3745727 
^ 24696000 + F ^ I ~ 210039480000 ~ 945 ~ 198450 

2^ I 54861581223623 4108 3745727 
+n/r^C^ I " 420078960000 + l35^' ~ ^96900"^^ ' ^^"^^^ 

^ ^ ^ , 37259291367883 19766 j_ 19766 . 1573589 . 

I pGrG/l I 1 D4 — C4 H C3 

' 38887309440000 945 105 6804 ^ 

200739467 \ 2/^3817101976847353531 39532 19766 

^ 28576800 j + ^ ^ 8166334982400000 945~ ^ ^ 105 

80980811, 497748102211 \ ( 740566685766263 632512, 



297675 ' 6001128000 / ' " \ 17013197880000 8505 

36241943 . ^ ^2 ^ f 4763338626853463 39532 



3572100 



^ I ^ P ^ \ 34026395760000 1215 

C2I (8.79) 



. 3952556872585211 4180 4180 21723277 
^ ^ ^ ' ~ 340263957600000 ^ 189" ~ ^T^^ ^ 85050 

559512437 \ _ /^ 1008729211999128667 8360^ 4180, 
71442000 i ^ ^ I 2041583745600000 189 ^ 21 

85539428, 131421660271, \ ^o^/^ 393938732805271 133760, 



297675 ' 1500282000 ' / ^ \ 8506598940000 1701 

19247947 \ f 2523586499054071 8360 

893025 ] + ^ [ 17013197880000 ^ 243 

19247947, , 
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10710275715721975271 48220 _ 48220 2873636069 
^ ^ ^ ' 452891327565600000 ^ 2079 ~ 231 ^ 10291050 

961673201 \ J. ^2 / 170291990048723954490137 _ 96440 ^ 48220 
^112266000^^ J + ^ I 328799103812625600000 2079 ^ 231 

10844970868 . 183261101886701 . \ « / 6080478350275977191 



36018675 ' 1996875342000 ' / \ 124545115080540000 

_ 1543040 2451995507 \ ^2^/38817494524177585991 
18711 ~ 108056025 J ^ ^ 1 249090230161080000 

96440^ 2151995507, \ 

H Ca C2 (8.81) 

2673 ^ 216112050 ^ j ^ ' 

/ 22309979286641292041 251264 251264 
^ ^ 603855103420800000 ^ 10395 ~ 1155 

283300123, 1210188619 \ ^^0/^177435748292579058982241 



935550 130977000 ' / " \ 328799103812625600000 

502528 251264 451739191 47705202493793^ 
" 10395 " ^ 1155 ~ 1440747 ^ 499218835500 

/ 6365809346912279423 8040448 512808781 
^ ^ ' 124545115080540000 ^ 93555 ~ 21611205 

T^C I 40517373495580091423 502528 _ 512808781 . 
^^•'^ ^ ^ ' 249090230161080000 ^ 13365 ~ 43222410 ' ^ ' 

126207343604156227942043 3387392 _ 3387392 
^ ^ ^ ' 2463815086971638400000 ^ 135135 ^ ~ 15015 

51577729507 2401246832561 
^ 158107950 ^ 243486243000 

^ 2 68296027149155250557867961293 _ 6774784 ^ 3387392 
^ ^ 122080805651901196900800000 135135 '^^ 15015 

79117185295, 108605787257580461, 
-Ca + .^^.^00^0. .^o.^^ C2 



243486243 ' 1096783781593500 



2^ [ 189306988923316881320303 108396544 
^ ^ ^ ' 3557133031815302940000 1216215 



90143221429 M / 1201733391177720469772303 



3652293645"^^ J ^ ^-^^^^^ I 7114266063630605880000 
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6774784 90143221429 \ 
^ 173745 ~ 7304587290 j ^ ' ' 

T C C 12032123246389873565503373 3498932 ^ 3498932 
^ ^ ^ I 181090408892415422400000 ^ 135135 ^ ~ 15015 

2288723461 106764723181157 
^ 6548850 ^ 10226422206000 

2 10076195142551036234891679659 6997864 3498932 
^ ^ 1^ 17440115093128742414400000 135135 15015 

81672622894 , 448416864235277759 



243486243 ^ 4387135126374000 

, 196243066652040382535303 111965824^ 93360116539, 
^ ^' 3557133031815302940000 1216215 ^^ 3652293645 ^^ 

I 1242840812874342588467303 6997864 
^ ^ ' 7114266063630605880000 ^ 173745 

93360116539 . 

' 7304587290 ' ^ ' ' 

. 994774587614536873023863 720484 ^ _ 720484 
^ ^ ^ ' ~ 12072693926161028160000 ^ 27027 ~ 3003 

6345068237 37428569944327 
^ 17027010 ^ 3408807402000 

2 72598193631729215117875463981 1440968 720484 
^ 122080805651901196900800000 27027 3003 



2101051892878 , 461388998135343407, 



6087156075 ' 4387135126374000 

^2^ J 40540032063650894708251 23055488 481761665447 
^ ' 711426606363060588000 243243 ~ 18261468225 

^ ( 256205552272074402170491 1440968 

^ ^ ' 1422853212726121176000 ^ 34749 

481761665447^^ ^ ^^^^^j 



36522936450 
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